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We study deterministic escape dynamics of the discrete Klein-Gordon model with a repulsive 
quartic on-site potential. Using a combination of analytical techniques, based on differential and 
algebraic inequalities and selected numerical illustrations, we first derive conditions for collapse of 
an initially excited single-site unit, for both the Hamiltonian and the linearly damped versions of the 
system and showcase different potential fates of the single-site excitation, such as the possibility to 
be "pulled back" from outside the well or to "drive over" the barrier some of its neighbors. Next, we 
study the evolution of a uniform (small) segment of the chain and, in turn, consider the conditions 
that support its escape and collapse of the chain. Finally, our path from one to the few and finally 
to the many excited sites is completed by a modulational stability analysis and the exploration 
of its connection to the escape process for plane wave initial data. This reveals the existence of 
three distinct regimes, namely modulational stability, modulational instability without escape and, 
finally, modulational instability accompanied by escape. These are corroborated by direct numerical 
simulations. In each of the above cases, the variations of the relevant model parameters enable a 
consideration of the interplay of discreteness and nonlinearity within the observed phenomenology. 
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I. INTRODUCTION 



a. The escape problem for systems of coupled objects. The escape problem for single particles of coupled degrees 
, of freedom or for small chains of coupled objects, concerns their exit from the domain of attraction of a locally 
^ ■ stable state [13, [HI- Its investigation is crucial for the understanding of evolutionary processes in natural sciences, 
\ representing the transition from a stable state to another, or to extinction and collapse. 

Transition state theory deals exactly with this kind of problem by referring to the case where the considered 
objects, initiating from the domain of attraction of a locally stable state, escape to a neighboring stable state crossing 
a separating energy barrier. A natural example of this sort involves a bistable potential, where the barrier corresponds 
to the height of the saddle point separating the two minima of the potential. This mechanism has been widely explored 



^1 , in numerous physical settings, including (but not limited to) phase transitions, chemical kinetics, pattern formation 
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and the kinematics of biological waves - see, e.g., Refs. [14|, |31L 123, |38l and references therein. Depending on the 
geometry of the underlying potential energy landscape which may possess multiple minima, topological excitations 
with a "kink shape" (which can be thought of as instantons) interpolating between the adjacent minima may exist 
| and play a critical role in the study of transition probabilities. 

■ Another example of such an escape process is the one whereby the considered objects, initiating within the basin of 
. . .' attraction of a locally stable state, escape towards infinity. This process is often termed collapse and it is associated 
with the emergence of finite time singularities for the corresponding dynamical system. Collapse or blow-up in finite 
time is quite widespread as an instability mechanism in diffusion processes and wave phenomena in nonlinear media 
0, HH, Ho-IHI ■ Generally, the collapse phenomenon features the domination of nonlinearity when the initial data are 
sufficiently large; however, the escape process leading to collapse may arise in more subtle cases where the initial data 
may be contained even deeply within the domain of attraction of the stable state. This process indicates the existence 
of interesting energy exchange and transfer mechanisms between the interacting objects. We remark that, generally, 
collapse does not occur in integrable systems related to solitons or intrinsic localized modes . As concerns discrete 
systems, the most prototypical among the known integrable differential-difference equations, namely the Ablowitz- 
Ladik and the Toda lattices, also do not exhibit collapse behavior [3, H3|. In both the continuum and the discrete 
cases this feature can be associated with the presence of an infinity of conservation laws, which seems highly restrictive 
towards the possibility of collapse-type events. 

A seminal work examining the underlying energy mechanisms supporting escape is the work of Kramers (23l . |30| , 
implying that external, mainly stochastic, energy sources generate energy dissipation and fluctuation enabling the 
escape process. These mechanisms have been investigated in many variants of damped and driven systems and the 
theory of thermally activated escape has been generalized to systems with many degrees of freedom, with numerous 
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physical applications [llj, [l8| . 

b. Deterministic escape for nonlinear Hamiltonian systems of coupled oscillators. The possibility of an unforced, 
purely deterministic Hamiltonian system realizing an escape process from the region of phase space sustaining oscilla- 
tions around a stable fixed point will be our focal theme herein. This nonlinear dynamics scenario has been investigated 
in Refs. [HI, Ojl, H(| for a system of coupled harmonic oscillators subjected to a nonlinear external potential, namely: 

ti„ - (£/„+i - 2U n + U n -i) + W'(U n ) = 0. (1.1) 

The system (|1.1|) . known as the Discrete Klein-Gordon (DKG) equation, is a fundamental model used in different 
contexts, including crystals and metamatcrials, ferroelectric and ferromagnetic domain walls, Josephson junctions, 
nonlinear optics, complex electromechanical devices, the dynamics of DNA, and so on 
In Refs. [lj| [20| the case of a cubic potential of the form: 

W{U) = ^U 2 -^U\ (3>0, (1.2) 

is considered. This is a prototypical example of potentials which may support collapse, since it possesses a metastablc 
equilibrium and becomes negatively unbounded after crossing its maximum. Importantly, it may support escape 
dynamics associated to collapse. Here, collapse will be taken to mean that one or more of the oscillators acquire an 
amplitude U — > oo. For instance, it is established [l9l[20| that an initially almost uniformly supplied energy, associated 
with an almost flat mode, may progressively be redistributed internally, creating unstable growing nonlinear modes, 
and eventually concentrated within a few units forming a localized mode that may overcome in amplitude the height 
of the potential well. More precisely, the whole chain is initialized by an appropriate randomization to an almost 
homogeneous state enabling non- vanishing interactions and the exchange of energy among the units. The total energy 
of this initial state is greater than the potential energy of the saddle which defines the depth of the well; however, the 
position of each of the oscillators within this state is close to the bottom of the respective well and thus the ener gy o f 
each unit is significantly below the energy barrier of the saddle. Due to the emergence of modulational instability [29(, 
localized excitations are formed and the total energy becomes concentrated in confined segments of the chain. This 
energy exchange mechanism results in the emergence of a critical mode in which one of the lattice units overcomes 
the saddle point barrier. 

c. Model and principal aims. In the present work, the collapse instability, as well as the above escape sce- 
nario will be examined by both analytical (based mainly on analysis techniques involving differential and algebraic 
inequalities for suitable norms of the system) methods and selected numerical experiments that will be used to com- 
plement/corroborate the analytical results. Our model will be the following Klein-Gordon chain, characterized by a 
quartic potential and a dissipative term: 

U n +jUn-(U n+ i-2U n + U n _ 1 )+u; 2 d (U n -pU^=0, (i > 0, 7 > 0. (1.3) 

The lattice is assumed to be infinite {n S Z) and we consider initial conditions: 

U n (0) and U n {0) € £ 2 , (1.4) 

thus the results refer to spatially localized solutions (as the definition of the phase space £ 2 directly implies). In 
some cases, especially in numerical simulations, we will consider Dirichlct or periodic boundary conditions. When the 
dissipation parameter 7 = 0, the DKG system Q1.3P describes the equations of motion derived by the Hamiltonian: 

-, +00 -, +00 2 +°° q 2 +°° 



while the case 7 > corresponds to the linearly damped analogue of the problem. Equation (JT73J) implies that each 
individual oscillator of unit mass evolves within the quartic on-site potential 

W (U) = ^-U 2 - (1.6) 

The repulsive ^-potential (|1.6p . illustrated in Fig.QJa), has a stable minimum at J7 m i n = 0, corresponding to the rest 
energy E min = W(0) = 0, and two unstable maxima located at C/jJ ax = T^j, corresponding to E m&x = W (^f-^=j = 

4/3- 
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System (|1.3p shares similar phenomenological properties with those of simplified models for DNA dynamics 0, [H, Q . 
To be more specific, the prototypical similarities are the following. The case j3 > is the physically relevant case 
for the study of energy-localization mechanisms in hydrogen-bonded crystals or DNA molecules, since in this type of 
models the potentials soften for large amplitudes [231 . In particular, the potential (|1.6p (which is also known as "soft 
quartic potential") may correspond to the hydrogen-bridge bond between nucleotides. As it is typical in the modeling 
of chemical bonds 0, Section 2.1, pg. 269], (|1.6[) represents qualitatively a potential with a hard repulsive part and a 
softer attractive part, and it diverges as U —¥ oo. The region [t/~ ax , U^ ax ], which corresponds to the potential well of 
depth -E max , describes the region of amplitudes below the breaking of the chemical bonds. 

Here, the escape problem can simply be described as escaping from the potential well (towards infinite values of 
the field) by crossing the saddle points U^ ax in the configuration space. In order to have escape of particles over 
the energy barrier E max to the region U > U max , a sufficient amount of energy has to be supplied. Regarding the 
large-time behavior of solutions, the escape leads to blow-up of solutions in finite time. 

d. Structure of the presentation and main findings. The paper organization and the main results can be described 
as follows. Our aim is to progressively examine the potential blow-up of initial conditions involving more oscillators. 

First, in Section II, we study the collapse problem for the Hamiltonian system (7 = 0), i.e., the formation of finite 
time singularities, and find the following. 

• Since the system can be seen to belong in the class of second-order evolution equations of divergent structure, it 
can be naturally treated, with respect to global non-existence, by the energy- type methods of (Ty| . In particular, 
using differential and algebraic inequalities we derive conditions under which suitable norms (such as the £ 2 norm) 
will diverge in finite time. 

• Focusing on initial conditions of a single excited unit, with the other units initially located at the minimum 
t^min = 0, we derive an analytical value for the initial position of the single unit outside the potential well, 
which serves as a sufficient condition (threshold) for global non-existence. This value depends strongly on the 
parameter uf^ controlling the interplay of discreteness and nonlinearity. Our numerical simulations demonstrate 
that, indeed, blow up occurs when the analytical prediction is satisfied and a concerted escape may follow 
when we are in a regime where the coupling effect is significant. This is the "drive over" phenomenon, where a 
particular oscillator drives its neighbors over the collapse barrier. A reverse phenomenon that we also identify 
numerically and explore herein is the "pull back" effect, where the attraction of an initially escaped unit by the 
neighbors back to the stability domain, leads to solutions which exist globally in time. 

• A detailed numerical analysis of the single excited site initial data reveals the existence of a "true" threshold 
for the position of the excited unit outside of the well, which acts as a separatrix for the above two different 
dynamical behaviors. For this threshold, the w^-dependent analytical prediction that we analytically identify (for 
collapse) serves as an upper bound; the relevant comparison of the analytical and true (numerically obtained) 
threshold is analyzed in the different limits of the discreteness parameter. 

• We also examine the global existence result of [l|| Lemma 2.1, pg. 455] in view of the Hamiltonian system 
()1.3D (7 = 0). We find the following: assuming negative initial time derivative of the £ 2 norm along with a 
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non-positive initial Hamiltonian energy fails to ensure global existence. The numerical results seem to fully 
support this concern. 

Next, Section III deals with the collapse problem for the linearly damped system (|1.3|) for 7 > 0. Using the abstract 
energy methods developed in Ref. HU (having the advantage that they take into account the geometry of the potential 
energy), we find the following. 

• We derive a 7-independent prediction for the position threshold outside the potential well for the single excited 
unit. The numerical experiments of section fill Al verify the relevance of this criterion and that, even for sig- 
nificantly large values of 7 > 0, the induced friction docs not modify the scenarios found in the conservative 
case. In particular, friction does not affect the "pull back" effect by increasing the threshold value (as might 
intuitively be expected). 

• The 7- independence of the threshold criterion, illustrates that only the strength of the binding forces are 
responsible for the "pull back" or "drive over" effects. To this end, the criterion is also applied successfully 
to the case where 7 = 0, yielding an alternative and improved upper bound to the threshold for collapse or 
stability. This is also corroborated by our numerical simulations. 

In Section IV, we continue our analysis by studying the escape dynamics for multi-site excitations of the Hamiltonian 
system (7 = 0). Here, we extend our arguments for a single-site excitation to the case of a few and ultimately to 
many site excitations, but this time, positioned inside the potential well. Our results in this setting are as follows. 

• We derive an analytical threshold for collapse based on the initial "energy" of a short-length lattice-segment, 
depending on the discreteness parameter u>j. The numerical results (based on a three excited units initial 
configuration) validate the theoretical expectations. The threshold separates the remaining of the segment in 
the potential well from its escape which, in turn, is leading to collapse. 

• Combining a violation of a small initial data global existence result (which can be proved by implementing the 
methods of Ref. @), with the modulational instability analysis (based on the discrete nonlinear Schrodinger 
(DNLS) approximation (29[), we reveal the existence of three different regimes for the amplitude of a plane 
wave initial configuration, distinguishing between the following dynamical behaviors: modulational stability, 
modulational instability without escape and finally to modulational instability combined with escape. We also 
examine how the violation of a small data global existence criterion connects to the numerical amplitude value 
beyond which modulational instability leads to escape. 

• For the linearly damped DKG chain, we perform numerical simulations which reveal a transient modulation 
instability regime for small values of damping. Nevertheless, the instability is completely suppressed eventually, 
in accordance with the predictions of the damped DNLS approximation. 

Finally, Section V offers a discussion and a summary of our results. The complementary appendices section contains 
the proofs of the global nonexistence and global existence results which have been used in this paper for the analysis 
of the escape dynamics. 

e. Notation. We shall use when convenient, the short-hand notation {U n (t)} n ^z = U(t) £ £ 2 , {U n (0)} n ^z = 
U(0) £ £ 2 and {U n (0)}n& = U(0) £ £ 2 , for the solution for t > and the initial conditions at t = respectively. In 
this notation, A2 stands for the one-dimensional discrete Laplacian 

{A 2 C/}„ eZ = U n+1 - 2U n + f/ n _i, 

defined on £ 2 . We shall also denote by 

+00 +00 
(U,Y) P = J2 U n Y n , \\U\\%= U l 

n— — 00 n— — 00 

the squared-^ 2 inner product and norm respectively. 

II. THE CASE OF THE HAMILTONIAN SYSTEM 

We may start discussing the conditions for the existence of finite-time singularities, and their relevance to the 
problem of escape dynamics. To this end, we have reviewed and implement in the discrete case the method of 
Ref. [l6|, on the derivation of a differential inequality for the norm x(t) := ||£/(t)||| 2 , and verify that it blows-up in 
finite time under appropriate sign conditions on the initial Hamiltonian and the time derivative of the norm. 



Theorem II. 1 We assume that the initial data for the system \1.3\) are such that H(Q) < and (U (0), U(0))jp > 
0. Then T max < oo. More precisely, the solution blows up (in the sense that \\U(t)\\p becomes unbounded) on the 
finite interval (0,Tb), with 

rp 9 l|t/(0)|| 2 , T , T 

Jb = 2 : , and J ma x < -ib- 2.1 

(U(0),U(0))p 

Proof: Sec Appendix A. □ 



A. Numerical Study 1: Connection of the results of Theorem HI7T1 with escape dynamics 

Initiating here our numerical studies, we intend to discuss the possible connections of Theorem lII.ll with the problem 
of escape dynamics. The numerical studies will consider the system (|1.3[) in a form involving the variable discretization 
parameter e > 0, namely: 

U n -e(U n+1 -2U n + U n - 1 )+co 2 d {U n -pU*) = 0, p > 0, e = ^ (2.2) 

on the interval [—L,L], supplemented with Dirichlct boundary conditions. The solution of (|2.2I) is described by the 
vector U G R K+2 

U(t) = (U (t),U 1 (t),...U K+1 (t)), U (t) = U K+1 (t)=0, (2.3) 

where U n (t) :— U(x n , t), x n = — ^ + nh, n = 0, . . . , K + 1, for all t > 0. With the change of variable t —> j^t, we 
rewrite the initial-boundary value problem for (|2.2[) in the form 



U n ~ (U n+ i - 2C/„ + tfn-i) + n 2 d (U n - /?C/3) = 0, t > n = 1,...,K, n 2 d = hW d , (2.4) 

U{0), U(0) e R if+2 , (2.5) 
[/o = E/jr+i = 0, t>0. (2.6) 

Our rescaling of the lattice parameter out of the problem makes it important that we add here a comment about the 
nature of the continuous and the so-called anti-continuous limit. In the latter, we need to consider within Eq. (|1.3|l the 
limit of uj 2 — ¥ oo. There, the nonlincarity dominates, and we are essentially in the regime of individual (uncoupled) 
oscillators. On the other hand, for small values of to 2 ,, the linear coupling term becomes important. However, due to 
the nature of our model, we are reaching the asymptotically linear limit (rather than the continuum limit) as u> 2 — > 0. 
We will perform numerical simulations for the simplest case of initial data 



U n (0) = AS nfi , (2.7) 

U n (0) - B6 n ,a, (2.8) 

where 5 n .o = 1 at site n = and zero elsewhere (Kronecker S), for A, B E M.. For the initial data (|2.7[) - (|2.8p . the 
condition for the initial Hamiltonian H(0) < reads: 

B2 , f 1 + ^l)a 2 < n 2 = hW d . (2.9) 



0, (2.10) 



2 V 2 / ~ 4 
Then, it can readily be observed that the quartic equation 

possesses the unique positive root (for A 2 ): 
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Therefore, for the initial data (|2.7j) - (|2.8p . the condition (|2.9[) is satisfied if 



A < -A» or A > A*. 



(2.12) 



It can be easily checked from (|2.11[) that 



-A, < XJ l 



1 



and A* > XJ l 




max 



VP 



max 



recalling that U^ ax denote the location of the saddle points of the on-site potential. The condition on the sign of the 
time derivative of the d 2 norm for the initial data (|2.7p - (|2.8[) reads 



First, we remark that condition (|2.12p is not directly related to the core question of escape dynamics, which is the 
escape from the potential- well of the on-site potential for initial configurations of the chain inside the well: assuming 
initially only one excited unit, the conditions of Theorem III.ll imply that the initial position of this unit should be 
located outside the interval (U~ ax , £/+ ax ) defining the well. 

On the other hand, exciting initially only one unit outside the well, with all the others located at the minimum 
U m i n = 0, is indeed of strong physical relevance in connection to another central question regarding the escape 
dynamics [13 . pg. 041110-4]: does this initially escaped unit continue its excursion beyond the saddle-point barrier or 
can it, in fact, be pulled back into the bound chain configuration U n £ (U~ a ^, t^max) by the restoring binding forces 
exerted by the neighbors? Alternatively, can the unit which is initially located outside the well, drag neighboring ones 
closer to or, in a more extreme scenario, over the barrier? The answer to these important questions should critically 
depend on the strength of the interactions between the oscillating units, imposed by the linear coupling. 

In view of the above questions, the conditions (|2.12l) - (|2.13p seem to be quite relevant, in the sense that they 
determine the behavior of the whole chain regarding the escape dynamics if the analytical threshold (|2. 12|) is crossed 
by at least one unit: in the anti- continuum limit e — > [or VL 2 d — > oo since, from our scaling Q 2 d = O (i) = 0(h 2 )], 
according to Theorem III.ll we expect that the initially escaped unit will continue its motion beyond the saddle- 
point barrier, while the other units should remain in the well. This is because in this case, the escaped unit is not 
interacting with the other units of the chain. The fact that ||f7(t)||#2 approaches infinity at finite time is due to the 
unboundedness of the energy of this unit in finite time; while the other units remain in the well, the excursion of the 
escaped unit implies actually the unboundedness of the sup-norm ||f7(t)||oo = sup n6Z |J7 n (i)| in finite time and, in 
turn, the unboundedness of the 1 1 C/ (t ) 1 1 £ 2 -energy due to the inequality: 



Note that in the infinite lattice only this side of the inequality holds in contrast with the finite lattice, where the 
equivalence of norms 



is valid in the ./V-dimensional space. 

In the discrete regime e = 0(1) (i.e., for moderate values of {I 2 , as e.g., D, d ~ 10 used in Ref. Q) and in the 
continuum limit e — > oo, it is expected from Theorem III.ll and (|2.14p that a concerted escape of the whole chain will 
take place, while the time of escape of the whole chain should depend on the coupling strength. 

The above analytical expectations and questions have been tested for fj = 1, h = 0.5 and varying the parameter 
ui d . The initial speed is taken to be B = 0.001, i.e., a small initial velocity of the one-excited unit is provided in order 
to satisfy the condition (|2.13p . 

The first observation of the numerical studies is the justification of \2. 13\) as sufficient conditions for blow-up 
and for the continuation of the excursion of the one-excited unit beyond the saddle point barrier when A> A*. 

The second question examined numerically, concerns the dynamics of the chain when the one-excited unit is located 
in the region 1 < A < A ifl i.e., beyond the location of the saddle point, but below the analytical critical value 
A*. There, our numerical simulations revealed that blow-up occurs for initial positions of the unit A < A*; hence, 
conditions (|2.12p - (|2.13p , although sufficient, are not strictly necessary for blow-up, which was found to occur for the 
one-excited unit even when its location is below the analytical value A* . 

However, and more importantly, the numerical studies revealed the existence of a new barrier (threshold) ^/thresh, 
located beyond the position of the saddle point [/jj ax = 1, with an important dynamical property. Although < 
^thresh, when the excited unit is initially located at position A satisfying [/jj ax < A < Wthrcsh < A*, it may be pulled 
back to the potential well and the solution exists globally. As mentioned above, the position of the barrier Wthresiu 



(U n (0),U n (0)) P = AB>0. 



(2.13) 



||tf(t)IU < \Mt)\\*- 



(2.14) 



TOH=o<||tf(i)||*< VN\\U(t)\\ 



(2.15) 




Figure 2: (a) Evolution of a single-excited unit initially located at A = 1.2, satisfying t/,i ax = 1 < A < A, = 2.45; the 
value of [/^ax is indicated by the dashed (red) line. Although the unit is initially beyond the saddle point (but below both 
the numerical threshold Wthrcsh ~ 1.8 and the analytical value A,), it is pulled back in the potential well, leading to globally 
existing dynamics. Lattice parameters are /3 = 1, fi^ = 1, for h — 0.5 and io\ — 4. (b) In this case, the single-excited unit, 
initially at A — 2, is pulling the adjacent neighbors towards escape from the potential well. Here, A, = 4.24 and Wthrcsh ~ 1-9. 
The lattice parameters are: /3 = 1, fij; = 0.25, for h = 0.5 and ui\ = 1. 



satisfying C/+ ax < Wthresh < ^4*, is determined by the restoring binding forces. This fact is demonstrated in Fig.OJa), 
showing the evolution of the chain in strong coupling regime (ft 2 , = 1). A single-excited unit, of initial amplitude 
A = 1.2 (satisfying t/+ ax = 1 < A < A*), is being pulled back inside the well, due to the sufficiently strong coupling 
with its neighbor units. This verifies our analytical prediction, that restoring forces (depending on the discreteness) 
can prevent the escape of elements of the chain initialized outside the potential well. This is what we will refer to 
as the "pull back" effect. After this effect takes place for the single-excited unit and it gets "retracted" inside the 
potential well (J7~ ax , E^m ax )> the wn °l e chain performs globally existing wave motions. 

The fact that blow-up occurs for A > A* = 2.44, as well as the result of Fig. Ufa), confirm the existence of the 
threshold Wthresh and motivate a further numerical investigation to determine its exact location. The corresponding 
numerical value is found to be t/thresh ~ 1-8 < A* for the near-continuum case (fl 2 , = 1). Concluding, we have shown 
that: 

1. blow-up occurs when A > Wthresh! 

2. a single-excited unit is pulled back to the potential well when /7+ ax < A < Wthresh! 

3. finally and in connection to our analytical result, A* provides an upper bound for the exact value of Wthresh G 
(C/+ ax , A*) and blow-up always takes place, as predicted for A > A*. 

The dependence of the analytical blow-up threshold prediction A* on the discreteness parameter ui 2 d (for fixed values 
of f3 = 1 and h = 0.5), is depicted by the solid (black) line in Fig. [3] The respective plot for the numerically obtained 
value, ^thresh is also depicted -in the same plot- by (red) circles. Comparing the two results, we conclude that A* 
tends to be more accurate as an upper bound for the "real" barrier [/thresh, for small to moderate values of uj 2 ,. On the 
other hand, a constant -but reasonable- discrepancy between the analytical and the numerical value is found in the 
intermediate regime between moderate and large values of uj 2 (the latter corresponding to the near anti-continuum 
regime of essentially uncoupled nonlinear sites). This discrepancy can be explained as follows: approaching the anti- 
continuum limit as uj 2 — > oo, it is naturally expected that the threshold value Wthrcsh — > t^max = since, due to the 
increasingly weaker interactions between the oscillating units, the system (|2.2I) is asymptotically reduced to a single 
repulsive Duffing equation for the single-excited unit -see Fig. [TJb). On the other hand, it follows from (|2 . 1 1 1) that 
A* — > -^| as — > oo. Thus, a constant difference between the threshold values is found in this regime, namely 

A* — Wthrcsh ~ (f° r large values of w^), which is observed in Fig. [3l 

In the asymptotically linear limit of uJ\ — > {Q 2 d — > 0), Eq. (|2.11[) predicts that A* — > oo. This prediction can 
also be justified, given the absence (in that limit) of the nonlinearity that is responsible for the blow-up. In the 
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Figure 3: Theoretical predictions for the blow-up amplitude for a single site excitation, as a function of ojJ, compared to the 
actual numerical value Wthrcsh- The solid (black) curve corresponds to the analytical upper bound A» [cf. Eq. ()2. 1 and 
the dashed (blue) curve to the analytical value A+,* [cf. Eq. (|3.12[) ]. Bullets (in red) denote the numerical values Wthresh- 
An asymptotic discrepancy appears for A, as lJ\ — s> oo while the accuracy of increases as the anti-continuum limit of 

uncoupled sites is approached- see discussion in the text. The parameter values are: h — 0.5, ft = 1. 



corresponding limit, we observe numerically that W t hrcsh — > oo, in accordance with the analytical prediction A* — > oo. 
The limits A* , U t hi CS h — > oo as w\ — s- are shown in Fig. [3J 

Furthermore, as the above limit is approached and due to the increasingly more significant role of the interactions 
between the units, if the excited amplitude has crossed the barrier Wtiucsh for blow-up, it will enforce the other units 
to escape from the potential well. This is "drive over" the barrier scenario, clearly demonstrated in Fig. EJb), which 
shows the evolution of the chain for lattice parameters Q 2 d = 0.25, i.e., to ? = 1 for h = 0.5. The single-excited unit 
with an initial amplitude i/thresh ^ 1-9 < A = 2 < A*, pulls the adjacent neighbors, enabling them to cross the saddle 
point towards escape from the potential well. 

The accuracy of the analytical prediction for the escape threshold will be improved in the next section dealing with 
the linearly damped version of (| 1 .3|) . The improvement, also depicted in Fig. |3j with the dashed (blue) line will be 
discussed there, in detail. 



B. Discussion on the sign-condition (|2.13|) . 

We conclude this section with a discussion of [la . Lemma 2.1, pg. 455] for the lattice dynamical system (|1.3[) - 
(|1.4p . According to the abstract results of Ref. [ly], keeping the condition "H(O) < and violating the sign condition 
([/(0), U(0))e2 > should lead to global existence. In terms of the initial data (|2.T[) - (|2.8|) . the result [HI Lemma 2.1, 
pg. 455] implies that the initially escaped unit will be pulled back and return to the potential well. For example, 
under the assumption 



x(0)=2(U(0),U(0)) i 2 < 0, 
we will reconsider the differential inequality (|5 . 18[) in the form 

x(t) . 3 



> 



[x(t)] 2 ~ 2x(t) ' 

Integrating (|2.17[) in the interval [0,t] for arbitrary t e [0,T], the interval of existence, we find that 



(2.16) 



(2.17) 
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Figure 4: Time evolution of the amplitude Uo, Hamiltonian H (denoted by H(t) in the figure), norm x and time derivative of 
the norm dx/dt for a single- excited unit with initial data ()2.7I) - (|2.8|) with A = 1.8 and B — —5.5; the lattice parameters are 
Q? d = 10 (oj^ = 40, h = 0.5), P = 1, corresponding to a highly nonlinear discrete regime. The top left panel shows the excursion 
to infinity for the one-excited unit; the top right panel depicts the negative (constant) Hamiltonian Hit) = H(0) ~ 22.7; the 
bottom left panel shows the change of the slope of the norm x(t) and its increase to infinity after time t*\ the bottom right 
panel depicts the change of sign of x(t) after time t*. 



1 3 /"* 1 



> - / -rvda, 



x(t) x(0) ~ 2 J x(s) 

which can be rewritten as 



1 1 3 /"* 1 



, < — -, — r - - / ~^-rds. (2.18) 

x(t) ~ x(0) 2 J x(s) V ' 

Then, the assumption x(0) = 2(U(0), U(0)) < 0, as well as the positivity of the norm function x(t) = \ \U(i)\\^ 2 > 
for all t € [0,T] in the interval of existence [0,T] should imply that < 0, i.e., 

±(t) < 0, for all t G [0,T]. (2.19) 
Integrating (|2.18p once more in [0,T], we see that 

x(t) < x(0), for all t G [0,T]. (2.20) 
Letting t — > oo in (|2.20p implies that T max = oo and 

limsup||£/(i)||2 2 < 1117(0)111,. (2.21) 

t— ^oo 

Hence, the solution C7(t) = {U n (i)} neZ should be defined in [0,oo) and be uniformly bounded. 

However, the validity of (|2.18|) for all t G [0, T] -the interval of existence- should be put under question: Although 

the continuity of x(t) = 2{U(t), U(t))p < in [0, T], the assumption ac(0) = 2(C/(0), U(0))e> < and ([2~T5|) guarantee 
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the existence of t\ > 0, such that x(t) < for all t G [0, ti], this assumption and (|2.18|) do not guarantee that t\ = T. 
Recall that the inequality f|5 . 14[) implies that x(t) > for all t G [0,T], hence the function x(t) is increasing in the 
interval [0, T]. Due to this fact, the existence of a t 2 > t\ > 0, such that x(t) > for all t G [t2,T], cannot be excluded. 

Our concerns on the arguments of fltl Lemma 2.1, pg. 455] have been tested numerically for the initial data 
(|2.7[) - (|2.8p . If the above arguments were valid, the condition H(0) < on the initial Hamiltonian implying (|2.12[) on 
the position A of the single-excited unit, together with (|2.16[) which for the initial data (|2.7|) - (|2.8j) reads as 

(U(0), U(0))jp = AB < 0, (2.22) 

should imply the pull-back of the initially escaped unit inside the potential well (J7~ ax , t/+ ax ). 

In Fig. [4] we present the results obtained, after numerically integrating Eq. (|1.3[) . with initial conditions of the form 
of Eqs. ([^7f |) -(|2T5 ]) . with A = 1.8 and B = -5.5, thus satisfying ([2~T2"j) . and the condition (gUJ) , sincc ( u (P)i U(Q))p = 
AB = —9.9 < 0. The lattice parameters are f2^ = 10 (lj^ = 40 and h = 0.5), /? = 1, while the initial Hamiltonian is 
negative 'H(O) = —22.7 as required, and remains negative and constant in the interval of existence, due to conservation 
of energy, as shown in top right panel of Fig. |4j In the bottom right panel of Fig. [4J we show the time evolution of 
the function x(t) = 2(U(t), U(t))p. The initial value of x(t) is negative, as required by [l(| Lemma 2.1, pg. 455], and 
remains negative up to some finite time, say t* . On the other hand, at t = t* (where x(t*) = 0), the function x(t) 
becomes positive, and remains positive {or all t > t* , justifying the contradiction with the arguments of Ref. fldl ] for 
global existence. Also, since x(t) > for all t £ [0, T], x(t) is strictly increasing, and the norm x(t) is concave up for all 
t G [0, T], as shown in bottom left panel of Fig. 0J Choosing a time t 2 > t* for which x(t) > for all t > t 2 > t* > 0, 
the results of Theorem III. 1 1 come into play again: since the system (|1.3|) - (|1.4|) is autonomous, using the time t 2 as an 
initial time, and (U{t 2 ), U(t 2 ))i2 as initial data, we may repeat the arguments of Theorem III. II establishing that the 
initially escaped unit can never return in the bound-chain configuration, as observed in the top left panel of Fig. |4j 

III. THE CASE OF THE LINEARLY DAMPED SYSTEM 

In this section, we will consider the conditions for collapse in the case of the linearly damped analogue of Eq. (|L .3[) , 
namely, 

U n + jUn-(U n+1 -2U n + U n _ 1 )+LJ 2 d (U n -0U^) = O, 7 >0, /3>0, (3.1) 
with the initial conditions: 

U n {0) and U n (0) G £ 2 . (3.2) 

In the damped system (|3.ip , we may consider the abstract methods of Ref. [35[ and prove global non-existence by 
replacing the sign condition on 'H(O) with an appropriate smallncss condition. Added to the assumption of the possibly 
positive Hamiltonian, a condition for a sufficiently large initial size of the quartic term of the Hamiltonian is derived, 
replacing the one -imposed in the Hamiltonian case- concerning the sign of the time derivative of the norm of the 
initial data. Of primary interest here, will be the discussion of the effect of the damping in the behavior of the whole 
chain, regarding its single node or possibly concerted escape from the potential well. 
The global non-existence result of the section is stated as follows. 

Theorem III.l We assume that the initial Hamiltonian satisfies 

n(0) < |j, (3.3) 

and that, initially, the quartic term is sufficiently large, in the sense that the i^-norm of the initial position satisfies: 

i 

\\u(o)\u* = ( E u ^°n ( 3 - 4 ) 



Then, the solution of LS. ip -i CO)) does not exist globally in time. 
Proof:. Sec Appendix B. □ 
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A. Numerical Study 2: Connection of the results of Theorem Ull.ll with escape dynamics 



As in section III Al in the numerical simulations we will consider the linearly damped version of (|2.2[) involving the 
linear coupling parameter e > 



U n + 'yUn-e(U n+1 -2U n + U n - 1 )+ul(U n -pU*) = 0, p > 0, 7 > 0, e = i 



(3.5) 



on the interval [—L,L], supplemented with Dirichlct boundary conditions. The initial-boundary value problem for 
the damped system (|3.5|) can be written with the change of variable t —> ht as 

tin + lhU n - (U n+1 - 2U n + U n -i) + n 2 d (U n - pUl) = 0, t > n = 1, . . . ,K, Q 2 d = h 2 u 2 , (3.6) 

17(0), fj(0) e R K+2 , (3.7) 

U = U K+1 =0, t>0. (3.8) 

In the numerical simulations we will again consider the single- unit initial excitation (|2.7[) - (|2.8[) . For the initial data 
(l2T7) - (pT8j) . conditions (|3~3]) - (pOJi are implemented as: 



Here, the quartic equation 



4/3' 



1 



£ 2 ft 2 



4 A { + 2 A 2^4/3 



= 0, 



(3.9) 
(3.10) 

(3.11) 



has in terms of A the two roots: 



A 2 ± , =2 



1 + % 

/30 2 



± 



1 -j — ^ 
1 ' 2 



Z? 2 



5 2/3fi 2 4/3 2 



Note that A\ „ > 0, + > „ and + > if and only if B € [~-^§p, • Furthermore, it can be seen 



(3.12) 



that 



A\ . lf > -, for all u)%,h,P > 0. 

Then, both conditions (|3.9p and (|3.10p are satisfied, as required, if 

A < -j4 +i * or A > 



(3.13) 



(3.14) 



As in the Hamiltonian case, conditions (|3.14j) imply that the initial position of the excited unit should be located 
outside the interval (C/~ ax , C/+ ax ) defining the well, due to (|3.13j) . We point out that this condition of blow-up bears 
no signature of the damping parameter 7 and, hence, can be used for the Hamiltonian case as well, a point evident 
in Fig. [3] and one to which we return below. 

However, the question of whether the single-unit, initially placed outside the saddle-point barrier, will escape or 
whether it will be pulled back into the well by its neighboring units, becomes even more interesting due to the presence 
of the damping: the question is if damping has an additional effect in increasing the value of the threshold between 
"pull back" and collapse. Since the threshold is increased as the strength of the coupling forces is increased (cf. 
Fig. |3l), one might expect that that the same would happen with the damping force. 

The time evolution of the damped chain, for lattice parameters SI 2 = 10 (us 2 = 40, h = 0.5 -corresponding to a 
moderately discrete regime) and /3 = 1, is shown in Fig. [SJ For this set of parameters the analytical threshold for 
the amplitude of the single-excited unit, calculated by p,12[) . is A + _* = 1.36. To test this prediction, we consider 
an initial amplitude of A = 1.4 and an initial speed B — as per the Theorem IIII.ll In order to investigate the 
friction- induced effects, we consider a relatively large value of the damping parameter, i.e., 7 = 2000. Figure [5] shows 
the escape of the excited unit to infinity. On the other hand, Fig. IHJa) visualizes the rapid increase of the kinetic 
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Figure 5: Evolution of the damped chain for lattice parameters Q? d = 10 (wjjj = 40, h = 0.5) and (9 = 1. The central excited unit 
n = has initial amplitude A — 1.4 and initial velocity B = 0. The large damping, characterized by the parameter 7 = 2000, 
does not affect the threshold value for the "pull back" effect. 



energy of the single unit n = 0, despite of the fact that the total energy is dissipated according to the identity (|5.25|) 
as long as the solution exists. The evolution of the kinetic energy is also shown in Fig. [6th), but for lattice parameters 
closer to the asymptotically linear regime, i.e., for Vt 2 d = 1 (and (3 = 1 as well). Comparing the two panels of Fig. HI 
it is evident that the initially excited unit at n = does not affect its neighboring units in the discrete case (a), in 
the sense that the energy remains localized at this site; on the other hand, in case (b), due to the stronger coupling, 
the central excited unit pulls the adjacent units at sites n = ±1 towards the saddle point, and their kinetic energy is 
increased compared to case (a). Numerical simulations have been also been performed for other values of 7, and they 
all have justified the prediction of Theorem IIII.ll when A > the collapse does not depend on the value of the 

damping parameter; furthermore, the true threshold Wti ir0 sh is also 7- independent. Therefore, the "drive over" and 
"pull back" effects are chiefly controlled by the strength of the coupling forces and not from the damping. The failure of 
the intuitive expectation that dissipation should decay away the motion of the oscillators and hence be less conducive 
to collapse is evident here. This is predominantly due to the principal role of dissipation in decreasing the overall 
energy of the oscillator chain; yet, the latter scenario can be achieved by large amplitudes due to the nonlinearity, 
hence the concerted effect of energy decrease due to dissipation and its achievement by increased (squared) amplitudes 
due to nonlinearity gives rise to the finite time blow-up analyzed above. 

Another interesting feature is the following. The analytical value provided by (|3.12[) . shows an increased 

accuracy -in the parameter regime ranging from the moderately to the highly discrete (and the anti-continuum)- as 
an estimate of the numerical threshold thresh, below which the single unit will be pulled back by the exerted forces 
of the neighbors, and above which the chain configuration collapses. This improvement is due to the fact that the 
proof of Theorem IIII. 1 1 takes into more detailed account the geometry of the potential energy W(U). Since the proof 
of Theorem IIII.ll is independent of the damping parameter, and is also valid in the Hamiltonian case, we return to 
Fig. |31 where the dashed (blue) curve depicts A + ^ as a function of (for fixed B = 0, h = 0.5 and (3 = 1). The 
increased accuracy of [dashed (blue) line] towards the anti-continuum regime is explained directly from (|3.12[) 
since A + _* — > t/,„ ax = as u> 2 — > 00, and the constant discrepancy for A* asymptotically vanishes. Also, in the 

asymptotically linear case of <J\ — > 0, A +jSf — > 00, thus it correctly predicts the asymptotic behavior oft/thresh in this 
limit. For moderate values of wj, a comparison of (|2.11l) and (|3.12l) indicates that A* may be preferable to 
in this intermediate regime. Nevertheless, ^4+,* accurately captures both asymptotic limits and in the intermediate 
regime both A* and are expected to be least accurate due to the relevance of the terms omitted in the derived 
(algebraic and differential) inequalities. 
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IV. SCENARIOS FOR ESCAPE FOR MULTI-SITE EXCITATIONS 



We now move progressively away from the single site excitation scenario and toward settings where the energy of 
the chain becomes concentrated within a (wider than one site, yet confined) region (l9j . Even if initially the energy 
is equally shared among all units, after a certain time, the dynamics may lead to its redistribution, so that at least 
one of the involved units concentrates sufficient energy, to overcome the barrier defined by the saddle point of the 
potential. In this section, we will first consider such an escape scenario, focusing on the evolution of a small lattice 
segment. More specifically, we will consider the evolution of a few-site excitation, by means of the energy methods 
developed in the previous sections. 

Then, we will study escape dynamics of the chain, induced by the modulation instability mechanism [Tsj, [2S|. In 
the latter case, the initial excitation will involve a plane wave extending throughout the chain. We shall review the 
conditions for modulation instability of plane waves for (|2.2[) and how this instability may be related to a potential 
escape of the modulationally unstable excitation. 

Finally, deriving conditions for the global existence of initial data of sufficiently small energy, we shall investigate 
how the possible violation of the global existence assumptions may be connected to the escape mechanism, and 
formulate conditions for escape dynamics. This approach will consider solely an initial condition in the form of a 
plane wave. 



A. Energy methods on a lattice segment and escape dynamics 



In view of the energy methods discussed in the previous sections and motivated by the first localization scenario 
discussed above, we focus on the study of the evolution of a fixed segment of the lattice. In particular, we seek 
appropriate conditions for the "initial energy" of this segment, which may lead to an escape process. 

We consider again system (|2.2[) on the interval [-L, L] supplemented with Dirichlct boundary conditions. Here, we 
will study the evolution of a lattice segment occupying the unit subinterval [—5,5] C [— L,L] together with the first 
neighbors adjacent to the points —5 and |. 

The number of oscillators located outside the piece of the chain of unit length I = [— |, |1 is 



(f-i)(A- + l) 



L 



(4.1) 



where \x] = min{n € Z | n > x}, x € R. Then the number of oscillators included in the unit interval / is 



m = K + 2 - 6. 



(4.2) 
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Wc distinguish between two different possibilities for the cndpoints of / -cf. cases A and B below. 

Case A. The simplest case is when the endpoints of / are occupied by oscillators. In this case, the lattice spacing 
satisfies 

(m + 1)^ = 1, (4.3) 

and the unit interval I consists of A = m + 2 oscillators, the m included in / together with the two endpoints. 

We denote by /' the interval which comprises / and the two adjacent ones. The interval /' consists ofA + 2 = m + 4 
oscillators. The length of /' is, due to (|4.3p . 

L > = l + 2h = (4.4) 
m + 1 

We also assume that the endpoints of I' occupied by the neighbors adjacent to the endpoints of I, are located at the 
sites k and k + A + 1. Then, we may decompose the solution configuration vector (|2.3[) as U = U 1 ^ 1 + U 1 , where 

U L V = (U ,U 1 ,...U k . 1 ,U k ,0,...,0,U k+x+1 ,...,U K+2 ), (4.5) 
U 1 = (0,...,0,U k+1 ,U k+2 ,...,U k+x ,0...,0). (4.6) 

The initial conditions C/ (0), U(0) £ K^" 1 " 2 are decomposed similarly. Since the decomposition is linear, it follows from 
(|2.2p that the elements U 1 ^ 1 and U 1 satisfy, respectively, the equations: 

Ul - eA d U,{ + uj 2 d U,{ - puj 2 d Ui U 2 =0, n = 0,...,K + 2, 

U^ 1 ~ eA d U^ + ujjU^ 1 ~ Ul = n = 0, . . . , K + 2. 

However, taking into account the form of [/ L \ 7 given in (|4.6[) , the equation for U 1 can be written as 

U,{ - eA d U^ + LJ 2 d U^ - pu) 2 d U* U 2 =0, n = k + 1, . . . , k + A, (4.7) 
Ul = U k+X+1 = 0. (4.8) 



Relabeling for convenience, the system (|4.7[) - (|4.8[) can be considered on the interval I' of the A + 2 oscillators j = 
0, . . . , A + 1 as 

Zj - eA d Zj + Lu 2 d Zj - fiiojZf = 0, t > j = 1, . . . , A, (4.9) 
Z(0),Z(0) e K A+2 , (4.10) 
Z Q = Z x+1 =0, t>0. (4.11) 

For brevity of notation, we indicate the initial data and the solution of (|4-9[) - (|4. as elements of the set 

S P = {ZeR x+2 : Z Q = Z x+1 =0}. 

The system (|4.9[) - (|4.11j) on the segment /' conserves the Hamiltonian on i.e., 

A+1 A+1 2 A+1 „ 2 ]=\+l 

3=0 3=0 3=0 3=0 

for all t € [0,T max ]. Let us recall that the following inequality holds 

A+1 A+1 

e^Z^+i-Ztf^^Y^Zl (4.13) 

3=0 3=0 



for all Z £ iSj/ , where 
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is the first eigenvalue of the discrete Dirichlet Laplacian on I' , 



-eA d $j = j = 0,..., A, 

$0 = *A+1 = 0. 



Case B. In the general case, the end-points of the unit interval / are not occupied by oscillators. In this case, the 
interval I' of the cut-off procedure has length 



L' 



Ah, when h > 0.5, 

(m + l)h, when h < 0.5, 



and the first eigenvalue of the discrete Dirichlet Laplacian is 



j ^shr (f) , when h > 0.5, 

2(m+l) 



when h < 0.5. 



With these preparations in hand, we shall start the investigations for conditions on the initial data Z(0), Z(0) of the 
segment, which may lead to escape dynamics. 

Proposition IV. 1 Assume that the initial data Z(0) G Sp of j^.ffp - f^Tnp satisfy 



|Z(0)|| 4 >*^±p. (4-15) 



Then, there exists < t\ < T maX; such that the solution Z(t) G Si> satisfies 



\\\Z(t)\\ A > for all t G [0,t x ]. (4.16) 



Proof: By using (|4. 13[) and (|4.14[) . we may observe that the Hamiltonian Hi>(t) satisfies 



^)>^±^(X>*| -^E^, forall te[0,U (4.17) 



Working as in Theorem lIII.ll we now consider the function 



*(P) == - ^V- (4-18) 



having the unique positive maximum 



( Ml + "I? t V 7 ^! + g 



$(p,) = *> at p, = v m V d . (4.19) 



Assuming that the initial data Z(0) G 5/' are chosen such that 



4> P* - 7=— , 

the continuity of the norm ||Z(i)||4 implies that there exists < t\ < T max , such that 



\Z(t)\U > p* = v//il t ? ^ ., for all t G [0,ti], 



which is (|4.16j) . o 
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Condition (|4.15[) is not establishing escape by itself, since it is certainly satisfied with t\ = T max , when the additional 
restriction on the initial Hamiltonian energy Hjr(Z(0)) on I', 

U I ,(0)<Hp*)= ifX \ + R 'f r , (4-20) 

holds. For instance, we may repeat the proof of Theorem IIII.ll for the excited segment on /' and show that the 
solution Z(t) of (|4.9p - (|4.11|) cannot exist globally in time. 

With the aim to avoid the restriction on the initial energy (|4.20[) and see if the condition (|4.15l) suffices for escape, 
we consider -in the set-up of (|4.9p - (|4.1ip - the invariant region introduced in [33[ (see also Ref. j47l|). For instance, we 
consider the set 

W(Z) := {Z e S P : I(Z) > and J{Z) < d} , 
where the functional I is given by 

A+l A+l A+l 

X(Z) = eY^iZ.+i-Z^+^Y.^-P^Yl 2 ^ 

J=0 3=0 3=0 

and d = inf J^(Z) denotes the infimum over all Z £ Sp of the functional 

A+l 2 3= X + 1 o 2 x + 1 

j=0 j=0 j=0 

The analysis of Ref. (33j for the continuum limit suggests that the set W(Z) is invariant under the flow associated to 
^4.9ty - fiJ.ll\ l and solutions are global in time. More precisely, if the initial data Z(0) G W{Z), i.e., 

I(Z(0)) > 0, (4.21) 

and H P (0) < d, (4.22) 

then the solution Z(t) of (|43|) - (|4TTj) satisfies 

Z{t) e W(t), for all t G [0,oo), (4.23) 

i.e., T m ax = oo- On the other hand, a solution which blows-up in finite time should initiate from initial data violating 
the condition (|4.2ip . It is interesting to observe that in the discrete case, an initial condition which violates {4.21\ ) 
satisfies j^,15[ ). 

Proposition IV. 2 // the initial data Z(Q) £ Sp satisfy 

2(Z{0)) < 0, (4.24) 

then they satisfy |^._/5| ). 

Proof: Condition (|4.24p reads as 

^ d J2zf(o) > £ f(Z j+1 (0) - Z 3 {0)f +^Ez|(o). 

3=0 3=0 j=0 

Applying (|4~T3)) and (pT2l?j) for Z(0), we have 

A+l A+l A+l 

e £(z J+1 (o) - z,(o)) 2 + Z H°) > (mi + <4) E ^ 0*1 

j=0 j=0 3=0 




Thus, we have derived that 
which implies that 

i.e., (|4.15p is satisfied, o 



fal\\z{Q)\\\ >{ lll +Jh\\z 



> V 7 ^! + 
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B. Numerical study 3: Evolution of the lattice segment 



Propositions IIV.1I and IIV.2I indicate that (|4.15|) may suffice for escape dynamics, taking into account that the 
connection with this behavior shall be established if fr4-15\ ) may describe initial configurations of the segment Z(0) 
inside the potential well. The numerical investigation of this issue is the purpose of the numerical study of this 
section. With the change of variable t — > j-t, we rewrite the system (|4.9I) - (|4.11I) as 

Zj - A d z 3 + n 2 z, - m\z] = o, t > o j = l, . . . , a, n 2 = h 2 ^ 2 , (4.25) 

Z(0),Z(0) e K A + 2 , (4.26) 
Z = Z x+1 =0, t>0. (4.27) 

We will consider a set of parameter values which will span the whole regime, from the asymptotically linear to the 
anti-continuum limit. For a lattice spacing of h = 0.5 we fall in case A, where the endpoints of the unit interval I 
are occupied by oscillators. In this case, the unit interval consists of A = 3 oscillators, with the m — 1 included in I 
[see also (|4.3j) ] together with the two endpoints. The length of the interval I', which is composed by I and the two 
adjacent neighbors, has length V = 2 and consists of A + 2=5 oscillators [see also (|4.4p ]. Since the number of points 
is A = 3, we will examine the evolution of the 3-point segment, by considering the simplest case of initial data, that is 
zero velocities Z(0) = (0, 0, 0), and initial positions of the form Z(0) = (A, A, A). For these initial data, the ^ 4 -norm 
is 

ii*(o)iu=fE^(°>) = (£^ (0) ) = ip/) =3U 

From (|4.14p . the eigenvalue ^i(I') — 2.34315, and the condition (|4.15[) results in a critical value for the displacements 
A, depending on u>d as follows: 



A>A clit := Vft T " * . (4.28) 

As will be justified below, the critical value A cr n is physically relevant for the escape dynamics. This is because in 
the discrete regime, it yields initial segment configurations inside the well and has the physically expected behavior 
in the limiting cases of the asymptotically linear and of the anticontinuous limit. 

In the left panel of Fig. [7] we show the time evolution of the 3-unit segment in the moderately discrete regime 
(Q 2 d = 10), with an initial condition, Z(0) = (A, A, A) where A = 0.79 > A cr ; t . As shown in this panel, the segment 
initially performs one oscillation, but subsequently reorganizes so that at time £ ~ 1.5 the central unit (j = 0) has 
already crossed the barrier ^thresh ~ 1.8 and it will lead the chain to collapse. In particular, it can never return to 
the potential well, according to Theorem lII.ll due to the autonomous nature of the system at hand. 

The energy transfer from the adjacent units of the segment to the central unit (resulting in growing amplitude 
oscillations of this unit), is visualized in the right panel of Fig. [7] The potential energy stored in the 3-unit segment 
(light colored area) is progressively localized within the central unit. The dark area between the patterns visualizes 
the passing from the bottom of the potential well U m - m = (with a maximal kinetic energy), as the units oscillate. 

The same escape dynamics have also been confirmed numerically, for the set of parameters in the strongly discrete 
regime (fl 2 = 100), with the initial condition A = 0.77 > ^4 C rit- The collapse time in this case is t ~ 6.4, which is larger 
compared to the respective one in the moderately discrete regime. This is expected due to the weaker interaction 
between the units, which results in a longer time interval during which the energy is redistributed, the more discrete 
the lattice becomes. Besides that, in the strongly discrete regime, the depth of the potential well is considerably 
increased, which has, as a result, a longer time needed for the excited units to climb over the saddle points and finally 
escape (cf. Fig.rjJ. 

The analytically obtained threshold A cr it, for which the segment Z(0) = (A, A, A) will exhibit escape dynamics, 
is plotted as a function of the coupling parameter ui d in Fig. [8] -see the continuous (black) line; for comparison, the 
numerically obtained threshold is also presented in this figure by (red) dots. The observation that lim w 3_>. A crit = oo 
as U! d — > reflects the fact that blow-up is not possible in the asymptotically linear regime. The numerical simulations 
verified the existence of a critical value uj d min , as we approach the asymptotically linear from the discrete regime, 
below which escape is prevented, and the energy stored in the 3-units segment is dispersed along the chain leading 
to global existence. For large values of u> d we observe that the numerical threshold converges slowly to the limit 
lim^^oo A cr it — 0.75. 
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Figure 7: Left panel: escape dynamics for the 3-units segment Z(0) — (A, A, A) with A — 0.79 > ^4 C rit; clearly, the central unit 
j = is escaping from the potential well. Right panel: contour plot showing the evolution of the potential energy of the chain 
in the case of an initially excited 3-units segment; the plot demonstrates the energy transfer from the adjacent units to the 
central unit and its localization therein. The parameters used are — 10 (ojJ = 40, h — 0.5) and /3 = 1. 
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Figure 8: The analytically obtained threshold A CI n, as a function of the coupling parameter to^ (continuous (black) curve), 
against the numerical value A csc ((red) dots) for the escape of the 3-units segment. Other parameters used are h = 0.5, P = 1. 



C. Modulational instability mechanism 

In this section, we will briefly review initially the modulational instability mechanism (l9| , with the subsequent aim 
to investigate its connection with escape dynamics of plane wave initial data [2^. This will complete our program 
of examining blow-up for the chain starting from one site excitations and progressively passing to few and finally to 
many site excitations. 

The modulation instability mechanism concerns the instability of a plane wave of the form 

U n = ^oe 1 ^"-"*) + c.c, (4.29) 

under small perturbations of different wave numbers. When this instability occurs, the exponential growth of the 
perturbations creates localized excitations. The question at hand here is if a sizeable growth of the perturbations may 
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create a critical localized excitation which, in turn, may lead to escape dynamics. Let us recall first that in the linear 
limit /3 = of (|2.2[) , the existence of plane waves (|4.29p is associated with the dispersion relation: 

c 2 = 4esin 2 f^ + uj 2 . (4.30) 



For q = and q = it, (|4.30j) gives uj 2 = w 2 , and ui 2 = uj 2 , + 4e respectively. Thus, the linear spectrum has a gap of size 
uidi and is bounded by the frequency w 2 lax = u d + 4e. 

/n the nonlinear case /3 > 0, we seek for solutions of (|2.2[) of the form 



C4W = ^ n (t)e- <Wit +c.c, 0„GC, (4.31) 

where <p n (t) denotes the varying envelope. Substitution of (|4.31|) into (|2.2[) implies that the varying envelope 4> n {t) 
satisfies the equation 

4> n - 2icj d <f n - uj<p n = e(<f> n+1 - 2<p n + - w 2 d (f> n + /3ojj (3\(f> n \ 2 <p n + c/) 3 n e 2iuldt ) + c.c. (4.32) 

Under the assumption of the slow- variation in time of 4> n (t) with respect to the main oscillation at frequency i.e., 

4> n << LOd4>m in the rotating- wave approximation, we may keep only the terms proportional to e u ° dt \ this way, (|4.32p 
results in the DNLS equation: 

2iu) d </>n + e(<?Wi - 2<p n + 4>„-i) + 3/3ujj\(f> n \ 2 (j) n = 0. (4.33) 
We proceed by seeking solutions of (|4.33[) of the form: 

<j> n {t) = 4> e w ", 9 n (t) = qn- Qt, O € R. (4.34) 

Inserting (|4.34|) in (|4.33[) . we obtain 

2u d Cj^e lqn + e0 o (V - 2e lqn + e^C"" 1 )) + 3^|^ |Voe i?n = 0, 
which leads to the following dispersion relation for the frequency uj: 

2uj d u = 4esin 2 (1) - 3/3w 2 |0 o | 2 . (4.35) 



In summary, the validity of the single-frequency approximation leading to the DNLS equation (|4.33[) for the varying 
envelope 4> n (t) of the solutions f|4.31[) of the DKG equation (|2.2p , assumes that the gap frequency aid is large compared 
to the other frequencies of the system, 

4e«wj}, ( 4 - 36 ) 

and the nonlinearity is weak, in the sense: 

Hi « 1. (4.37) 

Having the restrictions (|4.36|) and (|4.37|) in mind, one can study the modulational instability conditions of (|4.34[) by 
considering the ansatz 

K = (0o + Xn)e l{e " (4.38) 
where Xn and ifj n are small perturbations of the amplitude and phase, respectively, i.e., 

|Xn(*)| « 0o, and \Mt)\ « K(t)\- 
Then, substituting (|4.38[) in the DNLS (|4.33[) . it turns out that Xn and ijj n satisfy the system of linear equations 

2oJdXn + e [(Xn+i - Xn-i) sin q + 4> a (-ip n+ i - 2t/) n + ipn-i) cos q] = 0, (4.39) 

-2uJ d 4>o4n + e [{Xn+l - 2Xn + Xn-l) COS q - (/) (lpn+l - i>n-l) sin <?] + 6/3ujj<plxn = 0. (4.40) 

Assuming that the perturbations Xn and ip n have the form of plane waves, namely, 

Xn(t)= X ^ Qn - m \ (4.41) 
Mt) = iP e l{Qn - nt \ (4.42) 
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where xo an d V'o ar e constant amplitudes, while Q and Q denote wavenumber and frequency respectively, we can 
derive from (|4.39|) - (|4.40p the following system for xoi V'Oj 



-2i(u)dVL — e sin Q sin q) 



-4e^ sin 2 ( 2 ) cos q 



-4e sin' 



2 



cos q + 6(3uj^(j)Q 2i<j)o(uj c i^l — esinQsin q) 



Xo 



(4.43) 



For nontrivial solutions, we require the determinant of the matrix in (|4.43[) to be zero, which leads to the dispersion 
relation: 



(4.44) 



(wdf2 — e sin Q sin q) = e sin I — I cos q I 4e sin 1 — 1 cos q — 6f3uj d (j) 



From this, the condition for modulation instability, i.e., for f2 with a nonvanishing imaginary part can directly be 
obtained from (|4.44[) ; this condition is expressed (for cos(g) > 0; if this sign changes, so does the sign of the inequality 
below) in terms of the amplitude <fio and the wavenumber Q as: 



(4.45) 



60a>2#j > 4esin^ I -J cosg. 
In the case of the mode with q = 0, the condition for 4>q leading to modulational instability is: 



h > 



2esin 2 (g) 



»0,MI- 



Condition (|4.46|) should be combined with restrictions (|4.36|) and (|4.37|) which, in our case, are: 

1 » 

We also remark that in the case of modes with < q < % , since 



i>l and J\ >> 4e. 



(4.46) 



(4.47) 



Q 



6/3c^ - 4esin^ ( ^ ) cosg > 6/3wj$ - 4e, 



instability of a plane wave occurs if 



*8 > 



2 e 



3 puj 3(5h 2 uj 2 d 3/3fi 



(4- 



where we have used the scaling t — > jt as before. 



D. Small data global existence conditions and their violation 

The last scenario which will be considered is based on the violation of small data, global existence conditions, and 
the investigation of the possible connection of this violation with the generation of instabilities. For the derivation of 
the conditions for global existence, which shall guarantee non-escape dynamics, we shall consider a discrete variant 
of the method of Ref. Q , and use the energy functional: 

E{t) :=\\U{t)\\% + \\U{t)\\\= E ^+ E {U n+1 -U n ?+ul Ul (4.49) 

n— — oo n— — oo n— — oo 

where ||t/||? 2 denotes the "linear coupling energy" -norm 

||C/|||:= E (U n+ i-U n ) 2 +0J 2 d E Ul 



71— — OO 



n — — oo 
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We shall also consider the quantities 



E+ = 



which are real and positive if £7(0) satisfies 



1±. 1-£E{0) 1 



/8 



(4.50) 



£7(0) < 



2 K 1 + ^l)' 

Then, the small data, global existence result is stated in the following Theorem. 

Theorem IV. 1 We consider the energy function l{4.49\ ) and we assume that the initial data h3.2j) have energy 



(4.51) 



E(0) < min ^ 1, ■ 



(4.52) 



*( 1+ 4s) 

Then the solution of fl.3\) exists globally in time (T max = oo), and satisfies for all t € [0, oo) the energy bound 

E(t)<E- (4.53) 

Proof: See Appendix C. □ 

On the other hand, if the smallness condition (|4.51[) is violated, i.e., 



£7(0) > 



o 2 



2/3(1 + 4 



(4.54) 



the inequality (|5.86|) is valid for E{t) £ [0, oo) and the possibility of unbounded solutions cannot be excluded. 

With this observation at hand, we may investigate the dynamics of the chain for initial data satisfying (|4.54[) to 
see if this violation condition may give an analytical value of amplitudes leading to blow-up and escape. Let us also 
remark that 



2/3 1 



2<>- 



> 1 when fl 2 d > [3(1 + V2). 



(4.55) 



E. Numerical study 4: Modulational instability and violation of conditions for global existence 



We conclude this section with the last numerical study, investigating the modulational instability conditions (|4.4(il) - 
(|4.47p . and the violation of small data global existence (|4.54[) - (j4.55[) . We will consider lattice parameters in the 
strongly discrete regime, namely, h = 0.5 and uj 2 = 400 (f2^ = u! 2 h 2 = 100), while for the nonlinearity parameter we 
will use the value /3 = 0.1. We intend to examine (|4.48[) . i.e., the conditions for modulational instability and stability 
for amplitudes tfio > </>o,mi and 4>q < (ftp. mi, respectively. For our analysis we will focus on the region of modulation 
instability in the (q, Q)-planc defined in |29|, Fig. 2(b), pg. 3200], and particularly on the line q = 0; in this case, the 
corresponding domain for the wavenumbers is (0, and we fix Q = J. 

For the above set of parameters, the critical value of the amplitude </>o,mi defined by the right-hand side of (|4.48[) 



fo.Mi 



0.18. 



(4.56) 



The violation condition (|4.54[) will also be tested for plane wave initial data, and we will derive the relevant critical 
value of the amplitude $o,crit, as shown below. 

The plane wave (|4.29p can be rewritten as U n (t) = 2</>o cos(qn — u>t) and, thus, at t = 0: 



U n {0) = 20o cosqn, and U n (0) = 2uj<f>Q sinqn, 



(4.57) 
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In order to test the condition (|4.54[) numerically we should consider a finite lattice of K units. In this case, the initial 
data (|4.57[) have initial energy E(0) := E pw (0) given by 

K K K 

E pw (0) = 4lu 2 4>q sin 2 qn + 4<f>Q ^ [cos(g(n + 1)) — cosgn] 2 + 4(/> 2 ,£l 2 ! cos 2 qn. 

n—1 n—1 n—1 

It readily follows that the mode with q = has energy 

E pw {0) = 4Q 2 d 4>lK, (4.58) 
and the condition (|4.54[) will be satisfied for a chain of K units if 



00 > J- -pr ■= 0O,crit- (4.59) 

On the other hand, to get a A"-indcpcndcnt condition for the amplitude, we first observe that the mean value of 
the energy E pw (0) of the the q = mode denoted by E pw (0) is due to (|4.58|) 

1 



Since E pw (0) > E pw (0) the condition (|4.54[) will be satisfied for a chain of arbitrary size if 

^2 

£^,(0) > -y, 

2 ^( 1 + 4) 

implying the inequality 



£>o > 



i + 4 



*0,crit. (4.60) 



^ 8^ 

For the example of lattice parameters considered above, the critical value of the amplitude 'I'o.crit is found to be 

$ ,crit = l.H. (4.61) 

Note that both critical values for the amplitudes (|4.56j) and (|4.61j) arc physically relevant regarding escape dynamics 
since in the case (3 = 0.1 the saddle points arc located = T-^-g — 3.16 (and hence the initial excitation places 

all nodes inside their corresponding wells). It should be also recalled that the analytical condition (|4.60p rigorously 
does not ensure global existence or escape (blow-up), since it simply violates the global existence condition (|4.52[) of 
Theorcm lIV.ll for plane waves. The numerical simulations below will examine whether data satisfying the analytical 
condition ()4.60[) ] may lead to escape dynamics or not. 

The initial condition used in the numerical study is of the form [cf. Eq. (|4.57[) for q = 0]: 

U n (0) = 20 o + ecos(Qn), e « (4.62) 

where the amplitude of the small perturbation is e = 0.005. 

Figure demonstrates the transition from modulational stability to escape dynamics for plane wave initial data, 
by plotting the evolution of the ratio 

V(t)=m^^§-, (4.63) 

n «n(0) 



where the quantity 



U n = + l(U n+1 - U n f + i±Ul - (4.64) 



is the energy density of the system, and "H„(0) corresponds to the initial energy density, at t = 0. The figure 
demonstrates the transition from modulational stability (left) to modulational instability (middle) and finally to 
self-organized escape (right) upon varying the amplitude <po, as is explained in more detail below. 
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Figure 9: Evolution of the maximum energy density, normalized to the initial energy density, V(t) [cf. Eq. 1)4.63 Jl ], of the 
lattice chain. The initial condition, in all three cases, is: U n (0) = 2<^>o + ecos(Qn), with e = 0.005, Q = 5. The three panels 
correspond, respectively, to modulational stability for <fio = 0.07 (left panel), modulational instability but without collapse for 
4>o = 0.15 (middle panel), and escape dynamics for <f>o = 1.25 (right panel). The inset in the right panel shows the evolution of 
the maximum amplitude among the U n , clearly demonstrating the escape at t ~ 4.5. 



The left panel of Fig. [0J shows the evolution of (|4.63|) for 0o = 0.07 < 0o,mi and justifies (for this choice) the 
modulational stability of the plane wave; we observe that V(t) shows a small fluctuation around its initial value 
V(0) = 1, thus depicting the stability of the configuration over the entire time interval of the numerical integration 
of the system (1000 time units). Note that the same behavior was found for other values of the amplitude <fio, up to 
0o < 0.15. 

The situation changes drastically for (f>o = 0.15: beyond this value, modulational instability starts manifesting itself, 
as shown in the middle panel of Fig. [0J The onset of the instability is characterized by a significant increase of V(t) 
from its initial value V(0) = 1 (an order of magnitude larger than that observed in the previous case), which implies 
exponential growth of the excited mode. We mention that 0o = 0.15 < </>o.mi = 0.18; this deviation is expected due 
to the approximate description of the initial equation (|2.2|) by the DNLS model (|4.33|) . The long-time behavior of 
the system, in this modulationaly unstable regime, is shown in Fig. [T0J (the parameter values, as well as the initial 
condition, are the same with the ones in the middle panel of Fig. [3]). The contour plot of the energy density % n 
[cf. Eq. (|4.64|) ] reveals that up to t ~ 2000 a periodic pattern is formed, due to the modulational-instability-induced 
exponential growth of the excited mode. This is a so-called breather lattice (see e.g. [28| for a relevant discussion 
of such waveforms in integrablc models where they exist in analytically tractable form). However, after t — 2000, 
and due to the highly nontrivial dynamical interactions of such modes coming into play, at this stage, we observe 
the formation of moving breathers, each of which acquiring a small fraction of the total amount of energy These 
breathers interact with each other and exchange energy; during this process, it is observed that some "prevailing" 
breathers (the ones with the larger energy) seem to "absorb" the energy of the breathers that they interact with. This 
process continues for a long time leading to a gradual coarsening of the chain dynamics. Such processes have been 
debated at considerable length for Hamiltonian systems (with and without £ 2 norm conservation properties); a recent 
example of such a discussion containing an account of earlier work can be found in [22J. In our long time dynamics, 
one can clearly identify three strong "eventual" breathers; these are located far enough from each other, so that no 
interactions between them are observed (at least until the end of the simulation, at t = 10, 000). A similar effect has 
also been observed in Ref. Q (see Fig. 7 of this work and also references therein). 

The relevant evolution also brings forth the potential for connections of case examples such as that of Fig. HOI with 
chaotic dynamics. In particular, what we can see here, but also in Figs. 2 and 7 of pi, is that the MI manifestation 
leads not just to a simple unstable wavenumber, but rather to a band of unstable wavenumbers. The dominant one 
among them leads transiently to the formation of a structure resembling a breather lattice [28| , but as shown in both 
of the above examples, as well as in Figs, such as Fig. 3 of (28|, such states are unstable. Thus, the excitation of 
all the additional unstable modes eventually through an apparently chaotic evolution mixing these modes, leads to a 
self-organizing end result whereby a few dominant/robust discrete breather structures finally prevail. 

Let us now discuss the connection of the modulational instability mechanism with the emergence of escape dynamics. 
Generally, modulational instability gives rise to an exponential growth of the plane wave perturbations and, as 
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Figure 10: Contour plot, showing the long time evolution of the energy density "H n , for the same initial condition as in the 
middle panel of Fig. [5] Up to t ~ 2000, the chain follows a periodic (breather lattice) pattern due to the exponential growth of 
the initially excited mode. Beyond that point, we observe the creation of small breathers which collide and exchange energy. At 
the end of the simulation (t = 10,000), the total energy density is chiefly distributed among only 3 distinct localized breathers. 



a result, a large-amplitude localized mode -composed by one or few oscillators- may be formed. However, the 
modulational instability mechanism is not enough by itself to ensure escape dynamics, and modulationally unstable 
states corresponding to localized, even chaotic excitations, may exist globally in time. Our numerical simulations 
have shown that there exists a threshold of the amplitude 4>o separating the modulational instability regime from the 
escape dynamics regime; this threshold was found to be (for the set of parameters used in Fig. </»o = 0.98. The 
right panel of Fig. |H] depicts the evolution of ()4.63[) for cf>o = 1.25, i.e., beyond the threshold value. It is observed that 
for a short time interval V(t) remains close to its initial value while, at a later time, a sudden and sharp increase 
appears; this is also associated with a sharp increase of the maximum amplitude along the chain, leading to the escape 
from the potential well of the critical unit possessing this amplitude. The inset in the right panel of Fig. [5] shows the 
evolution of max„ |J7 n (t)| crossing the saddle point barrier. 

It is interesting to observe, that the analytical value 3>o,crit [cf. Eq. (|4.61[) ] appears to yield a value proximal to the 
amplitude threshold separating the modulational instability from the escape dynamics regime. We have confirmed 
numerically that an initial condition of the form of a randomly perturbed plane wave, namely, 

U n (0) = 2$ ,crit + e • noise, e = 1(T 4 , (4.65) 

leads to escape dynamics and blow-up in finite time, and the evolution of V(t) is similar to the one shown in the 
right panel of Fig. |H1 Additionally, it should be noted that our simulations have also confirmed the fact that $o,crit 
is almost independent of the number K of lattice sites, as seen in Table |TJ it is clearly observed that the numerically 
found critical amplitude takes the constant value 0o = 0.90 for chains composed by more than K = 100 sites (and up 
to K = 1000), while it takes a slightly larger value for smaller chains. 

The linearly damped system: Transient modulational instability within the exponential decay of solutions. In the 
case of the linearly damped system Q3.5P , seeking for solutions (|4.3ip under the slow- variation in time approximation, 
we derive the linearly damped analogue of the DNLS (|4.33[) , which is 

iuj d 4>n + iv d -4>n + e(4>n+i - 2(j) n + <j) n -i) + - /3uj}\<f) n \ 2 4> n = 0. (4.66) 







$o 


crit — 


= 1.11 






K 


20 


50 


100 


200 


400 


1000 


4>o 


0.98 


0.98 


0.90 


0.90 


0.90 


0.90 



Table I: The numerically found critical amplitude <j>o for escape dynamics for various numbers K of lattice sites. Here, 
"!?o,crit = 1.11 is the corresponding Jf-independent analytical prediction. 
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Multiplying (|4.66j) by <j> n and summing over n (in the infinite lattice or in the finite lattice subject to Dirichlet or 
periodic boundary conditions), we derive the equation 




implying the decay of the £ 2 -norm 

11^)11^=^*110(0)11^, (4.67) 

for all solutions of the approximating DNLS equation (|4.66[) . Furthermore, (|4.67[) and the norm relations ()2.14|) and 
(12 . 15[) imply for both the cases of the infinite or the finite lattice, the sup-norm decay estimate 

II^)IIL<^||0(O)||, 2 2 . (4.68) 

Hence lim^oo ||0(t)||oo — 0; an d all solutions of (|4.66[) decay at an exponential rate, independently of the initial data. 
However, the spatially uniform decay estimate (|4.68|) does not exclude a possible transient growth of perturbations, 
and consequently, a transient modulational instability behavior within the frame of exponential decay. This transient 
modulational instability is especially expected to be noticeable for small values of damping. 

The above suggestions have been tested numerically, by considering the evolution of the damped DKG chain (|3.5[) 
with perturbed plane wave initial data (|4.62[) . The values of parameters for h, j3, and Q, are the same as in the 
numerical study performed for the Hamiltonian system 7 = 0. The amplitude is 4>q = 0.15, inducing modulation 
instability in the Hamiltonian case. The two top panels of Fig. [IT] show the contour plots of the energy density H n 
of the dissipative system for 7 = 2 (left panel) and 7 = 0.02 (right panel), respectively. Although in both cases 
the solutions decay to zero as the approximating DNLS estimate (|4.68l) suggests, in the right panel we observe the 
formation of a transient breather lattice, which survives the exponential decay up to t ~ 100. The same breather 
lattice lives only up to t ~ 1, as shown in the top left panel, due to the stronger damping 7 = 2. The two bottom panels 
show the contour plot of the difference H n [U n \ 7> o] — H n [E/ n | 7 =o exp(— "ft)], i.e., of the energy density % n [C/ n | 7 >o] 
for the damped system (|3.5p from the energy density H n [C/ n | 7= o exp(— "ft)], where U n \ 1= $ denotes the solution of the 
Hamiltonian system (|2.2[) . In the bottom right panel the time evolution of the energy density difference is shown 
for the damping value 7 = 0.02. The peaks observed therein, within the time interval [20, 150], illustrate a transient 
modulation instability regime, where localization of energy and formation of a lattice of localized modes occurs. In 
the case of stronger damping 7 = 2, the transient instability regime is hardly observable. Numerical simulations for 
smaller values of 7 depict that the transient modulation instability regime is increased, as it is expected when the 
Hamiltonian limit 7 — > is approached. 



V. DISCUSSION AND CONCLUSIONS 



In this work, we have studied the escape problem in the discrete repulsive 4 -modcl. We have considered both the 
Hamiltonian (7 = 0) and the dissipative (7 > 0) variants of the model. We have assumed different types of initial 
conditions, ranging from single-excited units and small segments to plane waves. Our results address several relevant 
points related to the problem, as is briefly described below. 

In the Hamiltonian case, we have discussed how the system falls within the class of abstract second order evolution 
equations of divergent structure. Such systems can be treated, with respect to global non-existence, by the energy 
methods of Ref. [16j , involving the formulation of suitable differential inequalities enabling the identification of con- 
ditions for blow-up. For initial data in the form of a single-unit excitation (all other units being initially located in 
the minimum t/ m i n = 0), it is intuitively expected that escape occurs if the total energy of this unit is (sufficiently) 
greater than the energy barrier E maK , i.e., the initial position of the unit is located outside the potential well. Indeed, 
we have derived a suitable analytical value for the initial position of the single unit outside the well, which serves as 
a sufficient condition for the absence of global existence. However, these single-site initial data are very relevant for 
two important questions related to escape dynamics (already posed in Refs. (l9l.l2"oj): is it possible that the neighbors 
coupled to the escaping unit also get driven over the barrier and give rise to a concerted escape of the entire chain ? 
As an alternative scenario, can this escaping unit even be pulled back into the domain of attraction by the restoring 
coupling forces ? We have definitively shown that both the drive over and the pull back effects can be realized in 
numerical experiments of such chains. 

At this point, it should be remarked that the autonomy of the system allows to discuss the derived conditions in 
terms of the fate of the unit of the critical mode which, at a specific time, escapes from the well. In other words, 
time-invariance allows for posing as initial condition the position and momenta of the chain at an escaping snapshot, 
and consider as an initial time, any time referred to this barrier-crossing snapshot. 
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Figure 11: The top panels show the long time evolution of the energy density H n for the linear damped lattice 7 > 0. The 
bottom panels show the long time evolution of the difference H n [U n .\ 1> o] — H n \Un\-t=a exp(— ft)], i.e., of the energy density 
Hn [f/n|-y>o] for the damped lattice from the energy density Hn [t/„| 7= o exp(—yt)]. The left panels are for 7 = 2 and the right 
panels for 7 = 0.02. The initial condition is: U n (0) = 2cj>a + ecos(Qn), with e = 0.005, Q = ? and 4>o = 0.15. A transient 
modulational instability regime within the time interval [20, 150] is evident in the bottom right panel. 



The numerical simulations performed in Section [II Al justified that the analytical conditions capture the qualitative 
behavior of the chain and its dependence on the strength of the binding forces. Blow-up occurs when the analytical 
prediction is satisfied and a concerted escape follows when we are in the regime where the role of the coupling is 
significant. On the contrary, if the analytical condition is not satisfied then for suitable initial conditions, the initially 
out-of-wcll unit may be pulled back by the neighbors inside the stability domain, leading the initially escaping "spike- 
like" mode to give rise to a waveform which will exist globally in time. Furthermore, the numerical simulations 
revealed the existence of a "true" threshold for the position of the excited unit outside of the well, acting as a 
separatrix between the above two distinct dynamical behaviors. For this numerical threshold, the w^-dependent 
analytical prediction serves as an upper bound, with a justifiable error in the regime ranging from moderate to very 
strong discreteness. Both theory and numerics suggest that the threshold position moves far from the position of the 
saddles in the continuum limit. 

In the same setting, we have also examined the global existence result of [lfl Lemma 2.1, pg. 455], addressing the 
following question. Assuming negative time derivative of the £ 2 norm, while keeping the non-positivity condition of 
the initial Hamiltonian energy, does global existence arise for the solutions ? The answer we found is that the system 
under consideration [cf. Eqs. (|1.3I) - (|1.4[) ] serves as a counter-example for the aforementioned global existence result, 
a feature confirmed by our numerical simulations. 

We have also studied the linearly damped system [cf. Eq. (|1.3p with 7 > 0], for which the questions posed on the 
effect of the strength of the binding forces and their interplay with nonlinearity are additionally complicated by the 
role of damping. In that regard, we have tried to answer the following question: can the appearance of damping 
affect the position threshold which separates collapse from global existence ? To examine the above question, we have 
applied on (|1.3[) suitable modifications of the energy methods developed in Ref. [35| . These methods have enabled us to 
provide a 7- independent prediction for the position threshold of the single excited unit, and the numerical experiments 
of section 1111 Al have verified this independence. Theoretical and numerical evidence has also been provided for the 
improvement that this analytical estimate yields in connection to the true threshold. 

We have also examined the escape phenomenon for multi-site excitations inside the respective wells, trying to answer 
the following questions: is it possible to describe escape scenarios for such initial excitations ? Also, is it possible 
to derive quantitative conditions for these excitations (e.g., initial positions, initial amplitudes) which are sufficient 
for escape ? These questions were answered in the positive by using suitable modifications of our energy arguments. 
Finally, for plane wave excitations, relevant answers were provided by the analysis of the modulational instability 
mechanism. 

Motivated by the fact that escape is a phenomenon characterized by the concentration of energy within confined 
segments of the chain, in Section I IV Al we have examined the evolution of a lattice segment of unit length. With an 
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appropriate cut-off argument, we have derived a condition on the initial "energy" of the segment (in terms of the 
£ 4 -norm). Although this condition formally is not establishing collapse by itself, heuristic arguments based on the 
violation of the invariance principles of potential wells [33|, |4?| , indicate that it is relevant for initiating escape in the 
discrete regime. The numerical results validated the theoretical expectations and enabled us to recover the parametric 
regimes on which the escape process may be observed for our finite segment of excited initial data. 

Finally, we have examined escape, in terms of the instabilities of plane wave initial conditions (Section IV). In that 
regard, we have reviewed the derivation of the modulational instability condition for the amplitude of plane waves p9| . 
based on the DNLS approximation of the slowly- varying envelope solutions. On the other hand, as explained above, 
analytical energy arguments and direct numerical simulations revealed the existence of three different regimes for the 
plane wave amplitudes: modulational stability, modulational instability without escape, and modulational instability 
accompanied by escape. For the linearly damped DKG chain, the numerical results verified a transient modulation 
instability regime. The instabilities will be completely suppressed ultimately, in accordance with the predictions of 
the damped DNLS approximation. 

It is crucial to remark that modulational instability analysis structurally cannot be used for the examination of the 
self-organized escape dynamics for the DKG chain; the latter, is associated with blow-up in finite time (equivalent to 
the escape time when the initial data are in the potential well), while the modulational instability analysis relies on 
an approximation based on the use of the DNLS equation (|4.33[) ; for the latter, solutions always exist globally in time, 
independently of the initial data and the strength of the nonlinear term (26l . |4(| . This is a vast difference between the 
DNLS and DKG systems. Regarding the collapse behavior, the DKG demonstrates some analogies with its continuous 
counterpart, i.e., the nonlinear KG partial differential equation (PDE), which also may exhibit collapse depending on 
the size of the initial data and the strength of nonlinearity. Such analogies between DNLS and NLS equations are 
limited only to one spatial dimension and restrictively to the case of the cubic nonlinearity; in the case of focusing, 
power-type nonlinearities f{z) = \z\ 2p z, solutions of the one-dimensional NLS PDE may blow-up in finite time when 
p > 2, depending on the size and type of initial data Q, even in the linearly damped case [45j . Furthermore, in view of 
the modulation stability analysis, it should be recalled that the DNLS approximation drastically modifies the results 
concerning the stability domains, which may be deduced from a continuous NLS PDE approximation. In particular, 
it has been illustrated in (29[, that a stability analysis based on the continuous NLS approximation (the analogue 
of Eq. (|4.44|) for q, Q << 1) may fail, by erroneously predicting stability in regions of modulation instability which 
are correctly detected by the DNLS analo gue . Modulation instability effects appear in the NLS PDE equation with 
dissipation as it was shown, e.g., in Refs. [37l . |39| ; depending on the type and the size of the damping, modulation 
instability effects may be considerable. 

The results presented in this work may pave the way for future work in many interesting directions. Below, we 
briefly present a few such examples. 

• A natural possibility is to extend the present considerations to higher dimensional settings, where the role of 
the coupling will be more significant due to the geometry enforcing additional neighbors. 

• It is also relevant to consider via the methods and diagnostics presented herein the phenomenology of different 
types of potentials including e.g., the Morse potential, which is relevant to DNA denaturation as modeled by 
the Peyrard-Bishop model |3J] . We speculate that the techniques used herein may turn out to be more broadly 
relevant to problems involving hydrogen-bonds and, more generally, aspects of molecular dissociation. 

• While the present study has focused on cubic nonlinearities, it might also be of interest to consider quintic 
nonlinearities in the one-dimensional problem, or cubic ones in the two-dimensional case. In the present case, 
the DKG model is structurally closer to its continuum analog (regarding collapse properties) or to the DNLS 
(regarding MI features- but the DNLS model does not have collapse due to the conservation of the £ 2 norm). 



However, for the NLS equation in the continuum case, it is well-known 42j that a quintic nonlinearity is the 
threshold for collapsing dynamics in Id, while the cubic one is the corresponding threshold for 2d. Hence, for 
such models comparing/contrasting the collapse features of the lattice (DKG) model with the continuum (NLS) 
PDE problem would be of interest in its own right. 

• Finally, it is certainly relevant, following similar lines of approach as the work of fl9l . 20] to attempt to quantify 
the effects of noise in these systems and generalize in a probabilistic way the deterministic statements about 
escape presented herein for noise realizations with different correlation properties. 



Such studies are currently in progress and will be presented in future publications. 



28 



Acknowledgments 

J.C., B.S.R and A. A. acknowledge financial support from the MICINN project FIS2008-04848. The work of D.J.F. 
was partially supported by the Special Account for Research Grants of the University of Athens. P.G.K. gratefully 
acknowledges support from the U.S. National Science Foundation via grants NSF-DMS-0806762 and NSF-CMMI- 
1000337, from the U.S. Air Force via award FA9550-12-1-0332, as well as from the Alexander von Humboldt Founda- 
tion, the Alexander S. Onassis Public Benefit Foundation via grant RZG 003/2010-2011 and the Binational Science 
Foundation via grant 2010239. 



APPENDICES: Proofs of global non-existence and small data global existence results for the system (|1.3p 



In this complementary section, we provide the proofs of the analytical results concerning the global non-existence 
conditions, as well as the small data global existence conditions for the system (|1.3|) . 

We start our considerations from the local in time existence of solutions. Let us recall from [25l . pg. 468] that 
(|1.3[) - (|1.4[) can be formulated as an abstract evolution equation in £ 2 x £ 2 (see [1, [|| for the continuous analogue). For 

2 



instance, we may set A := A2 — uj 2 ,, and check that the operator 



B 



is a skew-adjoint operator, since 



B 



1 
A 



1 
A 



D(B) := 



D(B*) 



! 


X 







x,(e£ 2 , A x e 



(5.1) 



and D (B) = D(B*). The operator B, is the generator of an isometry group T(t) 
and bounded operators of £ 2 x £ 2 . Setting 



£(£ 2 x £ 2 ) -the space of linear 



U = 



u 
u 



, F(U) 



^P 3 



equation (jl.3| can be rewritten as 



U = BU + F(U) 



(5.2) 



In the above set-up, fixing T > and considering initial data U(0) = [U(0),U(0)] T G I 2 x £ 2 , a function U(i) 
[U(t),U(t)] T G G([0,T},£ 2 x I 2 ) is a solution of ([521), if and only if 



U(*) = T(*)U(0) + / T(t - a)F(U(a))ds := S(t)U(0). 
Jo 



(5.3) 



I 2 x I 2 and Eq. (|5.3[) can be treated exactly as in 



The nonlinear term defines a locally Lipschitz map F : £ 2 x t 
[26l . Theorem 2.1, pg. 94] in order to prove the following. 

Theorem V.l (a) For all U(0) G I 2 x I 2 , there exists T max (U(0)) > and a function U(t) G C([0, T max ), I 2 x I 2 ) 
which is for all < T < T max , the unique solution of ( f 5. 2\) in C([0,T],^ 2 x £ 2 ) (well posedness). 

(b) The following alternatives exist: (i) T max = 00, or (ii) T max < 00 and 

lim I |U(t) ||£2 X #2 = 00, (maximality) . 

(c) The unique solution depends continuously on the initial data: If {U„(0)}„ 6 n is a sequence in £ 2 x £ 2 such that 
U„(0) -> U(0) and if T < T max , then 5(t)U„(0) -> S , (i)U(0) in C([0,T],^ 2 x £ 2 ). 

Let us remark that for the local existence of solutions for p.ll) - (|3.2[) we have the same result of Theorem lV.il For its 
proof we just have to replace the matrix operator (|5.ip with the matrix operator: 



B, 



1 

A -7 



D(B) :- 



Zi lo€ £ 2 , Az G £ 2 



(5.4) 
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Appendix A: Proof of Theorem III. II 

We start by observing that the Hamiltonian can be written in the Lagrangian form 



H(t) = h\U\\%-V(U), (5.5) 



2 

where 

+oo 2 +°° a 2 +°° 



V{U):=-\ £ {U n+1 -U n f-^ Ul + ^L £ Ut (5.6) 



n— — oo n— — oo n— — oo 



The functional V : £ 2 -¥ R is differentiable (see 0, Lemma 2.3, pg. 121]), and for the derivative V'(U) at {/ G £ 2 , it 
holds that 



-t-OO + 00 +oo 



(V"(tf),Y>=- J2 (U n+ i-U n )(Y n+1 -Y n )-Lu 2 d J2 U n Y n + (3uj 2 d £ U 3 Y n , for all W £ I 2 . (5.7) 



71 — — OO 71 — — oo n — — oo 



Furthermore, by setting Y = U in (|5.7|) . we get that 

+oo +oo +oo 



(V'(U),U)v=- (U n+1 -U n f -uj 2 Y, UZ+fad E U n- (5-8) 



n— — oo n— — oo n— — oo 



With these preparations we proceed in various steps. 
Step 1: The inequality 

(V'(U),U)-W(U)>0, for all f/Sf 2 , (5.9) 
holds. Indeed, we see from (|5.6p and (|5.7p . that 

-i-oo +00 

(V' (17), U) P - W(U) = ]T C^f+i - ^n) 2 + "2 E ^ ^ °" ( 5 - 10 ) 



n— — oo n— — oo 



Thus the claim (|5.9[) is proved. 

Step 2: The function x(t) = \\U{t)\\ 2 2 satisfies the relations 

x(t) = 2(U(t),U(t))p, (5.11) 
x{t) > 2\\U\\l+8V(U). (5.12) 

Equation (|5.11[) follows by direct differentiation of x(t) = \\U(t)\\ 2 2 = J2n=-oo U n (t)U n (t). Note that all differentia- 
tions in the proof are justified by the local existence theorem of solutions in the interval [0,T max ) of local existence. 

For the inequality (|5.12[) we first differentiate (|5.11[) with respect to time and we substitute the expression U from 
the system (|1.3[) . which will be written for brevity as U = A2U — oj 2 U + f3u! d U 3 . Indeed, we have 

x(t) = 2(U,U) P +2(U,U) £ 2 

= 2\\U\\%+2(U,A 2 U-u; 2 d U + Pu 2 d U 3 ) i 2 

= 2\\U\\ 2 2 +2(V'(U),U) i 2. (5.13) 

Next, by using the inequality (|5.9[) of Step 1, which implies that (V'(U), U)p > 4V(U), we estimate the second term 
of the right (|5.13|) as 

x{t) > 2\\U\\%+2(V'(U),U) e 2 > 2\\U\\% +8V(U). 
Hence, the inequality (|5.12[) is proved. 

Step 3: Assume that the initial Hamiltonian is H(0) < 0. Then it holds that 

x(t) >6\\U\\% > 0. (5.14) 
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To prove (|5.14[) . we observe first that the equivalent expression (|5.5[) for the Hamiltonian implies that 

V{U) = \\\U\\%-U{t). (5-15) 

Due to the conservation of the Hamiltonian H(t) = H(0), (|5. 15|) can be rewritten as 

V(U) = l\\U\\%-H(0). (5.16) 

By replacing the term V(U) of the inequality (|5. 12|) proved in Step 2 with the right-hand side of (|5.16[) . we have that 

m>m\\i + %(\\\u\\i-nv))- (5.17) 



Since "H(O) < by assumption, from f|5 . 1 T[) we derive the inequality 

*(t)>2||J7||| 2 +8fi||{7||^ = 6||l>|| 2 2 > 0. 



Hence f|5.14[) is proved. 

Step 4-' Under the assumption H(0) < 0, the norm function x(t) = \\U(t)Wj 2 satisfies the differential inequality 

*(*) > (5-18) 

To prove (|5.18[) , we observe hrst that the Cauchy-Schwarz inequality implies the estimate 

[i(i)] 2 = [2(U, U) P } 2 < 4||t/|| f 2 2 ||t>||f 2 = ^{t)\\U\\l, 

which can be written as 

2 ^ [i(*) ia 



Then, from a combination of (|5.14[) and (|5.19p . the inequality (|5 . 1 8[) readily follows. 

Step 5 (finite-in time- singularity): In this final step for the proof of the theorem we shall use together with the 
assumption for negative initial Hamiltonian 'H(O) < 0, the additional assumption i(0) = 2(C/(0), U(0))i2 > (positive 
time derivative/increase of the £ norm initially). Under these assumptions, x(t) > in the interval [0,T max ) of 
existence and x(t) = ||E/(£)||f2 is strictly increasing in the interval of existence. 

Indeed, from the inequality (|5.14[) it follows that the function x{t) is increasing in the interval [0,T], T < T max . 
Since we have assumed i(0) = 2(E/(0), U(0) > and x(t) is increasing, it readily follows that x(t) > 0, on [0,T]. 
Therefore, x{t) is strictly increasing on [0,T]. The inequality (|5.18[) . can be rewritten as 

x(t)-2x(t)- (5 - 2U) 

Note that all the quantities involved are strictly positive on [0,T] as well as their initial values at t = 0. Integrating 
([5T2T)|) on [0, t] for t <T, we get that 

ln[x(t)] - ln[i(0)] > | (ln[x(t)] - ln[x(0)]) , 



which implies that 



x(t) > i(0) 



Integrating on [0, t] for t <T once more, we get that 



2 2 x(0) 

> -4nr*- (5-21) 



x{t)i x(0)i a:(0)i 
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Solving (|5.2ip in terms of x(t) we eventually derive the lower estimate 



^2x(0)f - x(0)i±(0)t 
The right-hand side diverges (blows-up) at finite time 

MO) n ll^(°)H a 



X(t) > ^ 2' ( 5 - 22 ) 



T h := 2- 



i(0) (17(0), 17(0))^' 

the time given in (|2.1|) . Therefore the left-hand side a:(i) = | [ C7 (t) | j 2 a cannot exist more than T\>. 

Appendix B: Proof of Theorem IIII. ll 

Step 1 (Energy identities): Multiplying (|3.1[) by U n and summing over Z (and by parts in the coupling term), we 
derive the energy identity 



1 d I 1 

2~dt 



{ 1 +oo 1 +oo 2 +°° o 2 +°° +°° 

5 E ^ + 5 E (^ + i-^) 2 + f E ^-^r E # +7 E *£ = °> ( 5 - 23 ) 
n— — oo n— — oo n— — oo n— — oo J n— — oo 

which, using the definition of the Hamiltonian (|1.5p . can be written as 

, +oo 

-H(t)+ 7 £ t/* = 0. (5.24) 

n— — oo 

Integrating in [0,t] for any i e [0,T], T < T max -the maximal interval of existence, we get the energy identity 

H(t)+j f \\U(s)\\%ds = H(0). (5.25) 
Jo 

Step 2 (Inequality for the Hamiltonian H(t)): For every t £ [0,T], the Hamiltonian satisfies 

i 

2 / +oo \ 2 2 +°° 



To derive (|5.26p , we first observe that 

E (^ + i-^) 2 + ^ # ( 5 - 27 ) 



1 +00 2 +oo 9 +oo 

n— — oo n— — oo n— — oo 



VJ1--00 / \n— — oo 



Furthermore, by applying the inequality 

i i 

/ +00 \ P / +00 \ q 

( E W) ^\ E W] ' forall 1 < 9 <p<c», (5.28) 

I — — OO / 

in the case p = 4 and </ = 2, we get that 



+ OC \ 4 / 4-00 



V n — — oo 



Thus, we have 

i 

/ +OO \ 2 +00 



E U n) < E Ul 
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With the above inequality at hand, we estimate the right-hand side of (|5 . 2T[) from below as 

+00 2 +00 2 +00 

\ e ( Un+1 -u n f+^ e m> ^ e u «- 

7i — — 00 n— — 00 7i— — 00 

. .2 / +°° \ 5 



> f E # ■ ( 5 - 3 °) 



\n- — 00 



The left-hand side of f|5.30p is the sum of the coupling and quadratic ?7 n -dependent terms of the Hamiltonian. Since 

1 



the kinetic energy term is always ^ J2n=-oo &n ^ Oj f° r au " * ^ [0; wc navc due to (|5.30p that 



«(*) = 5 E K + \ E (^ + i-^) 2 + f E ^ 2 -^ E 

n=— 00 7i=— 00 7i=— 00 

1 +00 2 +°° /Q 2 +°° 

>4 E (u n+l -u n f + "-± e ^ 2 -^ E ^ 



TL — — QO n — — OQ 71 — -— OO 

^ E E e 

\n— — 00 / n— — 00 

The claim (|5.26[) is proved. 

5*tep 5 (Derivation of the conditions i3.3\) and {3.$ and i/ieir consequences): The right-hand side of (|5.26p defines a 

functional for the scalar continuous function x(t) = \\U(t)\\ e i — fX]n=^-oo ^n) ■ This functional is of the repulsive 
</> 4 -type. For instance, (|5.26[) can be written in the short-hand notation 

H(t) > ^Wml - ^\\U(t)\\%, fo r all t G [0,T], (5.31) 
and motivated form the right-hand side of (|5.31j) . wc consider the function 

j( x ) : = ^ - ^A x \ x > 0. (5.32) 
Employing J, we will work out the conditions for x(t) = \\U(t)\\ei; we are interested on its unique positive maximum 

J{x m ) = g ati. = l (5.33) 
Now, let us assume that the initial data U (0) and U(0) are chosen to satisfy 

W(0) < J(s.) = ^ and ||i7(0)||,4 > a;, = -J=, (5.34) 

i.e., we assume that conditions p. 31) and (|3.4p are satisfied. Let us remark that the second of (|5.34j) and the continuity 
of x(t) = \ \U(t)\\ii, implies that there exists < t\ < T, such that 

\\\U(t)\\ £ 4 >x* = -^=, for all t G [0,ti]. (5.35) 

V P 

Due to the energy identity (|5.25|l . wc have the bound 

H(t) < H(0), for all t G [0,T]. (5.36) 
Then, from the first of (|5.34p and (|5.36[) . we derive that the Hamiltonian should be bounded by the maximum of J, 

uj 2 

H(t) < J(x*) = for all t G [0,T]. (5.37) 
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Then ([Qg]) . together with (jOTj) . imply that 

w 2 

J(x(t)) < H(t) < J(x.) = for all i e [0,T], (5.38) 

^P 

The consequences (|5.37p and (|5.38[) can be used as in [35], pg. 212], in order to prove that in (|5.35l) t\ = T, i.e., (|5.35[) 
holds for all t £ [0,T]. To see this, we assume (for the contradiction) that t\ < T. Then, there exists some t* > t\ for 
which 

*(**) = 111^(011* = *. = -^. 

Inserting this t* <G [0,T] in (|5.38p . the upper and lower bounds for H(t), 

uj 2 

J(a?.) = J(x(t*)) < H(t) < J{x*) = -J, for all t e [0,T] 
follow, meaning that J(x*) = %{t) for all i € [0,T]. This is a contradiction with (|5.37j) . Thus: 

|||tf(t)||<« > = -L, for all * € [0,7]. (5.39) 

VP 



It should be noted that the strict inequality of the assumptions p.3[) and (J3T4J) is crucial for the derivation of the validity 
of (|5.37|) and (|5.39[) . The consequence (|5.38[) has the geometric implementation in view of "Lyapunov functions", that 
under the assumptions (|5.34[) . both J and H calculated along the orbits {U, U} <G £ 2 x £ 2 are always less that the 
maximum J(x*) of the function J. 

Step 4 (Definition of a "blow-up" functional F(t) and the properties of its derivatives): We shall derive the differential 
inequality (|5.18j) for the function, 



F(t): = \\U(t)\\i 2+1 \\U( S )\\i 2 ds + 1 (T -t)\\U(0)\\i 2 +5(t + t o y 
Jo 

= (U(t),U(t)) e + 7 / (U(s),U(s)) i 2ds + 1 (To~t)\\U(0)\\ 2 2 + S(t + t ) 2 , (5.40) 



by applying the "quadratic form argument" of Pucci and Serrin [35|, cf. Theorem 1, pg. 206]. This argument will be 
explained in detail in Step 5. For its implementation it is necessary to discuss the derivatives F(t), F(t). In (|5.40j) . 
the constants to, To, 5 are positive and will be determined later. 

The functional F is well defined for all t £ [0,T]. The first derivative Fit) is found to be 

F{t) = 2([/,^ 2 + ||^ t ||[/( S )|| 2 2 d S |- 7 ||[/(0)|| 2 2 +2 ( 5(i + t ) 

= 2(U, U) P + ~f\\U(t)\\% ~ l\\U(0)\\% + m + to) 

= 2(17, U) P + y{(U(t),U(t)) P - (17(0), U(0))p} + 2S(t + t ) 

= 2{U, U) { 2+ / — {U{s) 7 U{s)) £ 2ds + 25{t + t ) 
Jo ds 

= 2(C/,Z7)^+ 7 y {{U(s),U(s)) i2 + (U(s),U(s)) l 2}ds + 25{t + t ). 
Therefore, the derivative F(t) is 

F(t) = 2(U, U)(2 + 2 7 / (U(s),U(s)) t 2ds + 2S(t + t ). (5.41) 



Differentiating (|5.41[) . we find a first expression for the second derivative F(t), 

F(t) = 2(17, tf )e- + 2(U, U) e - + 2 7 (17, U)<p + 26 



2 2 +2(U 7 U) e 2 +1 ^\\U\\ 2 2 +2S. (5.42) 



34 

Writing (|3.1[) in the shorthand notation U = A2II — jU — u) 2 d U + f3u>^U 3 , and substituting the derivative U in the 
second term of (|5.42p . we see that the expression for F(t) takes the form: 

F(t) = 2\\U\\% +2(U,A 2 U) P -2j(U,U) t 2 -2u 2 d (U,U) P + 2^(17, U 3 ) P + 7 ^||£/||f 2 + 2<5 

+ OO , + OO +OO . 



2||f>ll, 2 2 -2 2 (u n+1 -tr n ) a - 7 ^||tT||| 1 -2a;2 E t£ + £ U' n + 7^11^ + 2<J 

n— — 00 n— — 00 n— — 00 

+ OO +OQ +00 

-2 £ (U n+1 - U n f - 2u 2 Ul + *Po% E U n + 2S - 



n— — oo n— — oo n— — 00 



Hence: 

+OO -\-00 +OO 



On the other hand, from the Hamiltonian ()1.5j) . we see that the quartic term equals to 

+00 +00 +00 +00 

pu, 2 E ^" = 2 E tj n + 2 E (^+i-^) 2 + 2^ E U n-M(t)- (5-44) 

n=— 00 n— — 00 n— — 00 n=— oo 

Then, replacing the quartic term of (|5.43|) by the right-hand side of (|5.44[) , we write \F{t) as 

+00 +00 
-F{t) = ?,\\U\\l+ (U n+1 - U n ) 2 + cjj ]T U 2 -M(t) + 5. (5.45) 

n— — 00 n=— 00 

Recall that by the energy identity (|5.25l) . 

W(t)=K(0)-7 / ||t7(s)||f a ds. (5.46) 



By replacing the "H(£)-term of (|5.45|) with the right-hand side of (|5.46[) . \F(t) becomes: 

-F(t)=3||{7|||+ 2 (^ + i-^„) 2 +^ E U 2 -ilH(0)- 7 l \\U(s)\\%d S [+5. (5.47) 

n — — 00 n— — oo k J 

Furthermore, multiplying (|5.30p by 2, namely, 

i 

+00 +00 / +00 \ 2 

E (U n+1 - U n ) 2 + cjj E U n>"d( E ^ > forall ie[0,T], (5.48) 

n— — oo n— — oo \n— — oo / 

and from (|5.39[) . it follows that |||C/||| 4 satisfies 



/ +00 \ 2 2 

E ><^* 2 = ^ for all t€[0,T\. (5.49) 



Therefore, by estimating the right-hand side of f|5.4=8|) as it is suggested by ()5.49|) . we have 

-j-00 +00 



E (U n+1 - U n ) 2 + uj 2 UI>loW, for all t € [0,T\. (5.50) 



n— — 00 



Estimating further (from below) the sum J2n=-oo(^ n + 1 — U n ) 2 + lu 2 , J2n=-oo U 2 which appears in (|5.47[) by the right 
hand side of (|5.50[) . we eventually derive the following inequality for hF(t), 

±F(t)>3\\U\\j,+4v f \\U( S )\\%d S +{u 2 xl-4H(0)}+6. (5.51) 



35 



We require that the term uj^x 2 — 4FL(0) is positive. Hence, if 

u? d x\ - 4W(0) > 0, 
in terms of the initial Hamiltonian, the above requirement reads: 

«(°) < ^ = 2f = 

which is exactly the assumption (|3 . 3[) . Since this is satisfied as an assumption on the initial data, we may select 

26 = [C Ud xl - 4U(0)]. (5.52) 
With the choice (|5.52[) for 6, (|5.51[) can be rewritten as 



1 ■• • r* ■ 

=F(t)>3\\U\\%+4ry / \\U(s)\\j 2 ds + 36, 
* Jo 



2 

or equivalently 

F(t) > 6 (\\U\\ 2 £2 +6)+ 87 £ \\U(a)\\%ds. (5.53) 

Step 5 (Blow-up by non- continuation to [0, oo)): We argue by contradiction and assume that the solution U(t) of 
p.ip exists for all t € [0, oo), i.e., T max = oo. Then, at first, by choosing an arbitrary To € [0, oo) and since 6, to > 
(recall that a specific range for to will be fixed later), we have that 

F(t) > 0, for all t g [0,T ]. (5.54) 
By the expression for the derivative F(t) given in (|5.4ip . we see that at t = 0: 

F(0) = 2(U(0), U(0))p + 26t . (5.55) 
The choice of to will be determined by the necessity to assume F(0) > 0. This will be guaranteed by selecting 

to > -i([/(0),C/(0)), 2 . (5.56) 

Let us note that f|5 . 56[) is valid for any to > if (U(0), f/(0))^2 > 0, while the choice (|5.56l) for to is needed in the case 
where (U(0),U(0))p < 0. 

With the above choice for to, F(0) > and by (|5.53|) . which implies that F(t) > 0, we get that F(t) is strictly 
increasing in [0,Xb]. Together with the continuity of F(t), wc derive that F(t) > on [0,To]. Summarizing, we have 
derived that F and its derivatives are positive, i.e., 

F(t),F{t),F(t) > 0, for all t g [0,T ] for some T g [0,oo). (5.57) 

With f|5 . 53f) at hand, we may apply the "quadratic form argument" of Wh h recall that for any y) g R 2 , the 
quadratic form 

A X 2 +2B X y + Cy 2 , A, B, C <G R, A, C > 0, 
is positive definite (> 0), if and only if 

AC B 2 > 0. (5.58) 

The idea is to derive the differential inequality (|5.f 8[) for F(t) from (|5.58p . by proving that a quadratic form for 
suitable choices of real coefficients A, B, C involving F(i), F(i), F(t) is positive definite. 
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To apply this argument, we consider the quadratic form with the coefficients 

A := (U(t),U(t)) e 2 +7 / (U(s),U(s)) e 2ds + 6(t + to) 2 (5.59) 

Jo 

\%+! I \\U{ S )\\%ds + 5(t + t )\ 



Jo 

1 • • f* 

B := -F(t) = (U(t),U(t)) P+1 J (U(s),U(s))t*ds + S(t + t ), (5.60) 

C := (U(t),U(t)) t3 +7 / (U(s), U{s))fzds + S (5.61) 

Jo 

|2 / ||rV/„M|2 



= ITO + 7 / ||tf(*)||£ds + *. 

Jo 

Note that A,B > in [0,To], the first requirement for a positive definite quadratic form. Also, B > on [0, To]. 
With these coefficients, we have 

Ax 2 = (xU(t), X U(t)) P + 7 f\xU(s), X U{s)) P ds + X 2 S(t + i ) 2 , (5.62) 

Jo 

2B X y = F(t) X y = 2{ X U(t),yU(t)) e2 + 2 1 f ( X U( S ),yU(s)) P ds + 2 X yS(t + t ), (5.63) 

Jo 

Cy 2 = {yU(t),yU{t)) P + 1 f (yU(s),yU(s)) e 2ds + y 2 S. (5.64) 

Jo 

Before summing (|5.62|) - (|5.64[) to recover the quadratic form A X 2 + 2H X y + Cy 2 , we remark that 

0< \\ X U + yU\\% = ( X U + yU, X U + yU) P 

= ( X U, X U) P + ( X U, yU) P + (yU, X U) P + (yll, yU) B 

= ( X U, X U) e 2+2( X U,yU) e 2 + (yU,yU) e2 . (5.65) 

Summing now (|5.62[) - (|5.64[) , the right hand side of (|5.65[) appears as the sum of the inner product terms (outside and 
inside the integral). The sum of the remaining terms satisfies 

X 2 S(t + t ) 2 + +2 X yS(t + t Q ) +y 2 S = S ( X (t + to) + yf ■ 

Therefore, 

r t 

.2 , OTD,.„. , /-I..2 II TT I „. 7-7 1 1 2 , / \\„.TTf~\ I ..TTf~\\\2 



A X Z + 2B X y + Cy' = \\ X U + yU\\ z P + 7 / \\ x U(a) + yU(s)\fads + 5 ( x (t + t ) + y)~ > 0, 

Jo 

and therefore (|5.58|) should hold, i.e., 

AC > B 2 . (5.66) 

To derive the differential inequality (|5.18[) for F(t) with the help of (|5.66[) . we note the following: by comparing the 
definition of F(t) in (I5.40[) and the definition (|5.62[) of A we have: 

F(t) > A, for all t £ [0,T ]. (5.67) 

Next, to compare F(t) and C, we use (|5.53p and the definition of C (|5.64[) , showing that 

ft 

|2 , ex , o_, / llrV/„M|2 



F(t) > 6\\U\\j 2 +6S + 8j / \\U(s)\\j 2 ds 

Jo 

> 6\\U\\ 2 2 +65 + 6 7 / \\U{s)\\%ds = 

Thus, we have the differential inequality 



n 



^F{t) > C, for all t E [0,T ]. (5.68) 
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All the quantities in (|5.67p and (15.68[) are positive in [0,Tb], hence 

1 



F(t)F(t) > AC, for all t G [0,T ]. (5.69) 
By the definition of B in (f5JS0]> . 

B 2 = |[F(i)] 2 , for aU t G [0,T ]. (5.70) 

Since (f5T66|) and (|5l9]l . (|5J0| hold, we have that ±F(t)F(t) > j[F(t)] 2 , i.e., the differential inequality 

F(t)F(t) > ^[F(t)] 2 , for all t G [0,T ]. (5.71) 
The functions F(t),F(t),F(t) > in [0,T ] due to ([537) , and ([5T7T|) can be rewritten in the form of (|STTg|) . 

ff^ff. fo,all te[ 0,r ol . (5.72) 

Manipulating (|5.72l) exactly as in Step 5 of the proof of Theorem III. 11 we derive that the function F(t) cannot be 
continued further in time than [0, X 1 *), where T* is given by 

F(0) 

For instance, by choosing Tq = T* , the blow-up of the function F(t) in [0,Tb) contradicts the fact that To can be 
arbitrarily chosen in [0,oo) and that the solution U of (|3.1I) - (|3.2[) can be continued further than [0,Xb). 

Appendix C: Proof of Theorem IIV.ll 

Step 1 (Definition of the energy functional for the study of global existence): For convenience, we will denote the 
quartic term of the Hamiltonian (|1.5j) by 

a 2 +°° 

R(U(t)):=^ U n- ( 5 - 73 ) 

n — — oo 

Obviously, from the conservation of energy, we have 2H(t) = 2"H(0), and this can be rewritten by using (|4.49[) and 
(j5~73l) as: 

E(t) - 2R(U(t)) = E(0) - 2R(U(0)), for all t G [0,T max ]. (5.74) 



Step 2 (An inequality for R(U(t)) and E(t)): The energy functionals E{t) and R(U(t)) satisfy the inequality 

2R(U(t)) < -J^E 2 {t), for all t G [0,T max ]. (5.75) 
Zuj d 

To prove (|5.75[) . we first apply on (|5.73[) the inequality (|5.28[) for p — 4 and q — 2 to sec that 



n— — oo 



< ^W\\%. (5.76) 



On the other hand, the energy E(t) can be estimated from below as 



E(t)>uj J2 Ul (5.77) 
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4 2 



Recalling that \\U\\% = |(E^-oo U n) ' j = (E^-oo U n) > »* follows from^CS) that £(t) 2 satisfies 

#(*) > «5 ( £ 0*) =«5l 



which can be written equivalently as 



||C/||| 2 <^1. (5.78) 



Then, estimation of the ||?7||| 4 2-term in the right-hand side of (|5.76j) by (|5.78[) . results in 

^\ E 2 jt) _ p 



2R(U) < ^ ■ — ^ = -f^E 2 {t), for all t G [0,T max ], 



which is ([5775]) . 

iSiep 5 (Local in time estimates for E(t)—2R(U (i)) J: Some useful local in time estimates for the energy E(t) — 2R{U{t)) 
can be derived with the help of (|5.75[) . 

By writing the conservation of energy (|5.74[) as 

E(t) = E(0) - 2R(U(0)) + 2R(U(t)), for all t e [0,T max ], 

and estimating the 2 i?(C/ (t)) -term by using (|5.75|) . we derive a first estimate, namely 

E(t) < E(0) - 2R(U(0)) + J^E(t) 2 , for all t G [0, T max ]. (5.79) 



Note that (|5.75[) is valid for t = 0. Hence 

2i?(f/(0)) < ^£ 2 



Then, using again the conservation (|5.74l) and the positivity of R(U), we obtain 

E{t) - 2R{U{t)) = E(0) - 2i?(t/(0)) < E(0) - 2i?([/(0)) + 4R(Z7(0)) = E(0) + 2R{U(0)). 

Estimating the 2i?(£/(0))-term in the right-hand side by (|5.80[) . leads to the second -local in time- upper bound, in 
terms of E(0): 



E(t) - 2R{U(t)) < E(0) + -^£ 2 (0), for all t G [0,T raax ]. (5.81) 



Consequently, since f|5.81[) is valid at t = 0, the initial data satisfy 



E(0) - 2R(U(0)) < E(0) + -^E 2 (0). (5.82) 

2uj z 



Step 4 (Derivation of the condition on E(0) for global existence). We first assume that the initial data are small in 
the sense 

E(0) < 1. (5.83) 

Under condition (|5.83[) . we have that E 2 (0) < E(0), as well as that E 2 (0) < 1. Therefore, from the inequality (|5.82j) 
it follows that 



E(0) - 2R(U(0)) < E(0) + ^E 2 (0) < E(Q) + ^E(0), 



39 



or in other words, 



< E(0) ( 1 + JL ) . 



On the other hand, the estimate (|5.79[) can be written as 

/3 



< E(0) - 2R(U(0)) + i^E 2 (t) - E(t), for all t € [0,T max ]. (5.85) 



Estimating the E(0) — 2i?([/(0))-term in the right-hand side of (|5.85[) from above, by using (|5,84[> . yields the quadratic 
inequality for the derivation of conditions for small-data global existence 



©(£(*)) : = 7r^E 2 {t) - E(t)+E(0) 1 + > 0, for all t G [0,T max ]. (5.86) 



£%#w-m+m (1+4. 

The quadratic equation Q(E) = defined from (|5.86|) has the two positive roots (|4.50p . if -E'(O) satisfies the condition 
(|4"3T|) . Then, if (|4~5Tj) holds, inequality (|STM)) implies that 

E(t) G [0, £L] U +oo). (5.87) 

However, the function E(t) is continuous on [0,oo), therefore 

E(t) G [0, E-], (5.88) 

implying that solutions exist globally, and are uniformly bounded in time. 
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